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ABSTRACT 

The  sea  breeze  circulation  was  investigated  using  a  combination  of  acoustic  doppler 
sodar,  doppler  lidar  and  conventional  observations  in  central  California  on  the  coast  of 
Monterey  Bay  in  September  19S7.  The  study  was  called  LASBEX  (Land  Sea  Breeze 
Experiment)  and  used  the  combined  effort  and  resources  of  the  Naval  Postgraduate 
School  (NTS),  NOAA  Wave  Propagation  Lab  (WPL)  and  Naval  Environmental  Pre- 
diction Research  Facility  (XEPRF). 

A  monostatic  three-axis  phased-array  doppler  sodar  was  able  to  detect  the  various 
features  of  the  sea  breeze  front,  which  has  many  of  the  characteristics  of  a  gravity  head 
current  produced  in  laboratory  experiments,  such  as  the  longitudinal  vortex,  lobe  and 
cleft  structure,  and  strong  regions  cf  turbulence  and  convergence.  Profiles  of  the  wind 
field  anu  its  component  variances  were  generated  every  five  minutes.  Vertical  velocities 
up  to  1.0  m  s  were  detected  in  the  updraft  region  at  the  front.  Vertical  velocity  fields 
were  used  to  derive  the  horizontal  divergence  associated  with  the  sea  breeze  and  values 
of  convergence  in  the  range  of -4  x  10~3s-1  and  greater  were  calculated.  These  are  believed 
to  be  more  accurate  values  than  any  found  in  previous  sea  breeze  studies. 

Doppler  lidar  radial  velocity  information  was  useful  in  mapping  the  various  layers 
and  wind  structure  of  the  sea  breeze.  The  wind  field  of  the  sea  breeze  seen  by  the  lidar 
was  compared  with  the  sodar.  The  results  of  this  comparison  showed  that  the  sodar  was 
well  suited  to  investigate  features  under  600  meters  in  height,  such  as  the  updraft  region 
behind  the  front  and  the  variance  of  the  wind,  whereas  the  lidar  was  better  at  describing^ 
the  larger  scale  flow,  such  as  the  return  flow  of  the  sea  breeze,  the  convergence  zone 
along  the  sea  breeze  front  and  prevailing  synoptic  winds.  Together,  these  instruments 
gave  a  comprehensive  and  complete  description  of  the  sea  breeze  frontal  structure,  evo- 
lution and  movement. 
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I.      INTRODUCTION 

A.      GENERAL  DESCRIPTION 

The  land-sea  breeze  is  a  diurnal  mesoscale  circulation  that  occurs  on  sea  coasts  due 
to  the  differences  in  surface  temperatures  between  the  land  and  sea.  A  typical  horizontal 
temperature  gradient  is  1°C  20  km  (Atkinson,  1981).  This  situation  induces  a  thermally 
direct  circulation  that  alternates  between  onshore  or  sea  breeze  and  offshore  or  land 
breeze.  The  circulation  occurs  during  the  day  due  to  solar  heating  of  the  land  which 
warms  faster  that  the  sea.  This  results  in  a  horizontal  pressure  gradient  with  slightly 
higher  pressure  over  the  land.  The  pressure  differential  causes  air  to  flow  at  higher  levels 
towards  the  sea,  thereby  leading  to  convergence  aloft.  The  resulting  counter-flow  at  the 
surface  is  the  sea  breeze.  Fig.  1.  shows  a  lake  breeze  circulation,  which  has  a  similar 
structure  as  the  sea  breeze  circulation. 

The  land  breeze  is  due  to  the  same  thermal  structure  except  it  is  reversed.  At  night, 
the  land  loses  heat  quicker  than  the  sea  and,  therefore,  sets  up  a  circulation  opposite  to 
the  sea  breeze.  In  both  cases,  the  flow  at  the  surface  is  from  the  cooler  region  toward 
the  warmer  region.  By  comparison,  the  sea  breeze  is  the  stronger  circulation,  since  the 
heating  of  the  land  provides  for  a  more  intense  heat  source  to  drive  the  flow. 

The  sea  breeze  occurs  more  frequently  gi\en  a  weak  horizontal  pressure  field 
(Estoque.  1962:  Simpson  et  al,  1977).  This  is  because  the  thermal  effects  have  to  alter 
any  existing  conditions.  If  the  existing  pressure  field  is  strong  enough,  the  resulting 
gradient  winds  can  overwelm  the  land-sea  breeze  effect. 

Under  some  conditions,  an  offshore  gradient  wind  will  generate  a  frontal  zone  at 
the  leading  edge  of  a  sea  breeze  (Atkinson,  1981).  Estoque  (1962)  notes  in  his  model 
that  an  offshore  flow  counter  to  the  sea  breeze  will  enhance  the  horizontal  temperature 
gradient.  The  sea  breeze  front  is  also  characterized  by  gradients  in  humidity,  wind  speed 
and  wind  direction.  This  sea  breeze  front  is  similar  in  form  to  a  synoptic  scale  cold  front, 
because  the  cooler  moist  marine  air  is  advancing  into  the  warmer  dry  land  regime. 
Passage  of  such  a  front  is  a  measurable  feature  and  can  induce  significant  mesoscale 
meteorological  phenomena  such  as  showers,  thunderstorms  and  gusty  winds  (Pielke. 
1974;  John,  19S0;  Richiardone  and  Pearson,  1983).  The  horizontal  extent  of  the  sea 
breeze  is  large  enough  that  it  can  develop  a  component  parallel  to  the  coast  due  to 
coriolis  deflection  (Defant.  1951;  Walsh,  1974). 
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Fig.   1.      Typical  Lake  Breeze  Circulation,  (from  Keen  and  Lyons,  1978). 


B.      MOTIVATION  FOR  THE  STUDY  OF  LAND-SEA  BREEZE 

The  land-sea  breeze  circulation  and  sea  breeze  front  can  aflect  the  dispersion  of 
aerosols  and  pollutants  over  a  wide  area.  This  effect  has  been  observed  and  modelled 
many  times  (Lyons,  1972;  Lyons  and  Olsson,  1973;  Anthes,  197S;  Kondo  and  Gambo. 
1979;  Kitada  et  al.,  19S6;  Ogawa  et  al.,  19S6).  In  a  sea  breeze,  particulates  and  gases 
can  be  carried  inland  over  long  distances.  Collis  and  Russell  (1976)  summarized  some 
of  the  important  effects  that  particle  pollutants  have  on  climate,  such  as  cloud  forma- 
tion, precipitation  and  transfer  of  solar  radition.  Pollutants  can  enhance  cloud  forma- 
tion and  result  in  acidic  rainfall.  Particulates  and  gases  can  block  incoming  solar 
radiation  and  trap  outgoing  long-wave  radiation. 

In  the  case  of  lidar  (light  detection  and  ranging)  sensors,  concentrations  of 
pollutants  such  as  .\02  and  S02  can  attenuate  or  distort  the  return  signal  so  that  infor- 
mation may  be  misinterpreted  (Measures,  19S4). 


Measuring  the  intensity  and  mapping  the  structure  of  the  sea  bree/e  can  give  a 
better  understanding  of  how  particles  are  transported  and  dispersed.  Satellite  data  could 
be  used  to  detect  a  sea  breeze  occurrence.  In  addition,  forecasts  of  pollution  events  may 
be  feasible.  Surface  data  also  can  be  used  to  evaluate  satellite  observations  of  the 
planetary  boundary  layer  and  the  mesoscale  features  that  occur  during  the  sea  bree/e. 

C.      PREVIOUS  STUDIES 

A  study  was  conducted  by  Johnson  and  O'Brien  (1973)  of  the  sea  breeze  along  the 
Oregon  coast.  In  their  study,  the  sea  breeze  was  found  to  have  a  maximum  depth  of  700 
meters,  with  a  maximum  speed  of  7.5  m/s  at  the  surface.  The  zonal  wind  was  utilized 
to  describe  the  variability  of  the  sea  breeze  because  of  the  north-south  orientation  of  the 
Oregon  coast.  The  inland  penetration  of  the  sea  breeze  front  was  up  to  60  km.  A  return 
flow  was  found  between  1000  and  1500  meters.  Johnson  and  O'Brien  made  the  following 
observations:  the  low-level  onshore  flow  of  the  land-sea  circulation  was  located  entirely 
within  the  marine  layer;  the  marine  layer  tends  to  deepen  briefly  at  the  onset  of  sea 
breeze  due  to  the  low  level  convergence  ahead  of  the  sea  breeze  front;  the  return  flow 
above  the  inversion  appears  to  be  in  the  form  of  surges,  probably  in  response  to  the 
surges  observed  in  the  sea  breeze;  and  a  distinct  wind  maximum  follows  the  sea  breeze 
front  inland. 

Simpson,  et  al.  (1977)  observed  sea-breeze  fronts  that  penetrated  over  45  km  inland 
with  some  fronts  moving  as  far  as  SO  km.  These  fronts  had  an  average  advance  speed 
of  3  m  s.  They  also  showed  that  the  structure  of  the  sea  breeze  front  is  similar  in 
structure  to  gravity  currents  produced  in  laboratory  water-tanks. 

Mitsumoto  et  al.  (1983)  conducted  a  laboratory  experiment  which  used  a  water-tank 
to  simulate  the  sea  bree/e.  This  experiment  reproduced  many  of  the  characteristics  of 
the  sea  breeze  flow,  such  as  strong  vertical  motion  at  the  leading  edge  of  the  sea  breeze, 
compensating  flow  aloft  and  production  of  longitudinal  convective  vortices  within  the 
sea  breeze  flow. 

Ogawa  et  al.  (1986)  observed  the  gravity  head  structure  reproduced  by  Mitsumoto's 
laboratory  experiment  during  a  lake  breeze.  Ogawa's  experiment  utilized  a  kytoon- 
mounted  ultrasonic  anemometer-thermometer  to  measure  instantaneous  wind  and  tem- 
perature in  order  to  observe  the  turbulent  structure  of  the  lake  breeze.  He  found  that 
the  penetration  of  a  lake  breeze  produced  an  upward  rolling  motion  within  the  head  of 
the  front.  He  recorded  upward  flow  velocities  around  3  m/s.  Ogawa  found  that  the 
temperature  profile  had  peaks  indicating  the  frontal  passage.  These  peaks  were  due  to 


the  boundary  between  the  offshore  wind  aloft  and  the  lake  air  below.  Another  feature 
that  was  observed  was  the  development  of  the  Thermal  Internal  Boundary  Layer 
(TIBL).  The  TIBL  develops  due  to  the  temperature  discontinuity  between  the  land  and 
the  water  (Stunder  and  Sethurman,  1985).  The  TIBL  plays  a  significant  role  in  the 
fumigation  of  coastal  pollution  that  is  adverted  into  it. 

Estoque's  (1962)  model  had  a  maximum  penetration  of  sea-breeze  front  of  32  km 
and  a  minimum  penetration  of  18  km.  The  maximum  penetration  occurred  during  calm 
initial  conditions  and  the  minimum  penetration  occurred  during  offshore  counter-Cow. 

Pearson  et  al.  (1983)  used  a  model  and  found  that  with  a  mean  flow  perpendicular 
to  the  coast,  the  speed  of  the  sea  breeze  front  is  a  linear  function  of  the  onshore  com- 
ponent of  the  mean  flow.  Lie  found  that  the  maximum  vertical  velocity  in  the  front  at 
any  given  time  was  not  affected  by  the  mean  flow.  Another  characteristic  that  Pearson 
found  was  that  in  a  negative  shear  case  (upper  flow  is  from  land  to  sea),  there  was  an 
increase  in  the  vertical  extent  of  the  sea-breeze  front. 

For  a  general  descriptive  summary  of  sea  breeze  studies,  Atkinson  (1981)  provides 
a  comprehensive  review  from  many  locations  and  experiments. 

D.   LAND-SEA  BRE1  ZE  EXPERIMENT  (LASBEX) 

A  measurement  program  was  conducted  in  the  Monterey  Bay  area  in  September 
1987,  to  understand  the  structure  of  the  land-sea  breeze  effect  and  generate  data  to 
evaluate  many  of  the  models  which  have  been  used  to  demonstrate  the  behavior  of  the 
air  flow.  This  experiment  was  unique  because  it  utilized  a  varied  array  of  instruments 
including  doppler  acoustic  sounder  (sodar),  a  doppler  lidar,  rawinsonde  launches  over 
the  the  sea.  radiosonde  launches  over  land,  conventional  surface  recordings  of  temper- 
ature, humidity  and  wind  and  corresponding  observations  by  the  XOAA  9,10  and  GOES 
satellites.  This  combined  instrument  analysis  allowed  for  comparison  of  the  data  among 
many  sources.  The  final  result  will  be  a  comprehensive  view  of  the  land- sea  breeze  in- 
tensity and  structure.  This  study  focuses  on  data  from  one  sodar  supplemented  by  lidar 
observations. 

The  products  obtained  from  the  sodar  data  were  the  height  of  the  marine  boundary 
layer  and  time-height  cross  sections  of  wind  direction,  wind  speed,  vertical  velocity  and 
the  standard  deviation  of  the  vertical  velocity.  From  these  it  was  possible  to  compute 
time-height  cross  sections  of  horizontal  divergence  and  turbulent  kinetic  energy.  The 
time-height  cross   sections  will  be  compared  with  conventional  data   obtained  from 


radiosonde,  rawinsonde  and  surface  observations.    The  boundary  layer  height  observa- 
tions will  be  compared  with  the  time-height  and  conventional  data. 

The  lidar  provided  pictures  of  the  radial  wind  field  in  the  sea  breeze  front.  The  lidar 
observations  will  be  compared  with  the  sodar  to  find  differences  between  the  instruments 
ability  to  map  the  sea  bree/e  front  and  measure  the  associated  wind  field. 

The  combination  of  lidar  and  sodar  has  been  applied  before  to  analyzing  the  at- 
mosphere. However,  the  two  instruments  have  not  been  used  to  study  the  sea  breeze 
circulation.  Russell  et  al.  (1974)  used  the  combination  to  observe  the  atmospheric 
boundary  layer  in  an  urban  environment.  One  comparison  in  Russell's  study  was  the 
ability  of  either  instrument  to  detect  subsidence.  In  the  case  of  the  lidar,  detection  was 
degraded  because  aerosols  tend  to  be  excluded  from  descending  air.  For  the  sodar.  the 
temperature  gradient  associated  with  a  large-scale  subsidence  inversion  is  easily  detect- 
able. Lopez  (1977)  used  the  combination  to  study  the  properties  of  convective  activity 
and  small  cumulus  formation.  This  combination  is  effective  since  either  instrument  can 
detect  different  features  to  give  a  better  overall  picture  of  the  atmosphere. 

^  The  objective  of  this  thesis  is  to  evaluate  the  effectiveness  of  the  sodar  in  the  de- 
tection and  mapping  of  sea  breeze  events  and  to  compare  its  measured  wind  field  with 
that  of  the  lidar.  A  case  study  on  16  September  19S7  is  presented.  An  overall  look  at 
typical  characteristics  of  the  sea  breeze  in  the  Monterey  Bay  Salinas  Valley  will  be  gen- 
erated for  use  in  other  studies,  such  as  aerosol  transport  or  the  evaluation  of  satellite 
imagery. 


II.      DESCRIPTION  OF  EXPERIMENT 

A.  LOCATION  OF  THE  EXPERIMENT 

i  v 

LASBEX  was  conducted  from  between  15  and  30  September  1987  on  the  Monterey 

Bay  coast.    The  collection  of  data  started  at  1200  Pacific  Standard  Time  (PST)  on  15 

September  and  ended  at  1300  PST  on  30  September.    All  instruments  were  operated 

continuously  except  the  contracted  research  vessel  Silver  Prince,  which  was  stationed  in 

Monterey  Bay  to  monitor  surface  wind,  temperature,  humidity  and  pressure  as  well  as 

launch  rawinsondes.    The  ship  collected  data  from  1000  PST  until  1600  PST  on  various 

days  throughout  the  experiment.     The  instruments  evaluated  in  this  paper  were  the 

Pennylvania  State  University  (PSU)  sodar  and  the  NOAA  Wave  Propagtion  Laboratory 

(WPL)  lidar.   The  sodar  was  located  southeast  of  Castroville  at  the  mouth  of  the  Salinas 

Valley.     The  lidar  was  located  on  a  flat  knoll  approximately  1.5  km  from  the  coast, 

southeast  of  Moss  Landing.  It  was  on  a  line  of  310°  with  the  sodar,  which  approximated 

the  orientation  of  the  Salinas  Valley.    The  separation  between  the  sodar  and  the  lidar 

was  5.5  km.  The  Silver  Prince  was  stationed  in  Monterey  Bay  approximately  16  km  from 

the  lidar  location.   The  three  main  instrument  locations  are  shown  in  Fig.  2. 

The  terrain  around  the  sodar  site  was  at  elevations  ofless  than  30  meters,  with  ele- 
vations greater  than  30  meters  northeast,  starting  2.5  km  away.  The  interference  caused 
by  the  town  of  Castroville  was  considered  minimal  for  the  experiment.  The  terrain  sur- 
rounding the  lidar  was  generally  flat  with  higher  elevations  to  the  north  and  east.  To  the 
east  the  terrain  rose  to  levels  of  100  meters  and  greater,  up  to  400  meters.  To  the  north, 
the  terrain  generally  had  a  maximum  elevation  of  150  meters.  North-northeast  of  the 
lidar.  Elkhorn  Slough,  a  brackish  inlet,  interruped  the  elevated  terrain. 

At  the  lidar.  sodar  and  ship,  twenty-second  averages  of  vector  wind  velocity,  air 
temperature,  relative  humidity  and  surface  pressure  were  collected  at  the  stations  during 
the  experiment.  Radiosondes  were  launched  at  the  sodar  station  at  the  discretion  of  the 
research  team,  coinciding  with  significant  changes  in  the  environment. 

Satellite  AVHRR  imagery  during  LASBEX  was  available  from  the  Scripps  Institu- 
tion of  Oceanography. 

B.  ACOUSTIC  DOPPLER  SODAR 

During  LASBLX,  the  doppler  sodar  was  used  to  monitor  the  acoustic  backscattered 
energy  of  turbulent  fluctuations  in  the  atmosphere.    The  sodar  is  a  3-axis  monostatic 
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Fig.  2.  Location  of  Silver  Prince,  Lidar  and  Sodar:  The  Silver  Prince  was  at  lo- 
cation A.  The  Lidar  was  at  location  B.  The  Sodar  was  at  location  C.  Ele- 
vations are  in  feet.  Distances  are  in  kilometers. 


dopplcr  system  manufactured  by  Xontech  and  consists  of  a  phased  array  of  25 
transducers  which 'operated  at  a  frequency  of  1600  hertz.  The  pulse  repetition  interval 
is  10  seconds  per  beam,  giving  a  cycle  time  of  30  seconds  between  the  3  acoustic  beams. 
Two  cycles  were  integrated  per  minute.  A  wind  profile  was  sampled  every  five  minutes, 
giving  an  effective  number  of  10  cycles  per  record. 

The  sodar  uses  sound  energy  to  detect  the  small-scale  temperature  fluctuations  and 
discontinuities  in  the  mean  temperature  profile  (Crease  et.al.,  1977;  Koracin  and 
Berkowicz,  1988).  The  addition  of  doppler  capabilities  allows  for  the  recording  of  wind 
speed  and  direction  (Kaimal  and  Haugen.  1977).  The  backscattered  energy  provides  in- 
formation about  the  height  of  the  marine  boundary  layer,  wind  speed,  wind  direction, 
vertical  velocity  and  vector  winds.  N'eff  and  King  (1988)  demonstrated  that  the  sodar 
can  give  indirect  spatial  and  temporal  information  on  both  the  location  and  thickness 
of  turbulent  la\ers.  In  addition,  Neff  was  able  to  identify  mixing  processes  in  the 
boundary  layer.  Small  mixing  processes  could  be  attributed  to  turbulence,  where  as 
large  processes  could  be  associated  with  effects  such  as  Kcvin-Helmholtz  instabilities. 

C.      DOPPLER  LIDAR 

The  lidar  used  during  LASBCX  was  the  WPL  10.6  micron  doppler  system.  The  lidar 
operated  under  computer  control  and  had  real-time  data  processing  and  color  display 
capabilites.  The  system  could  compute  the  zero,  first  and  second  moment  of  the  echo 
spectrum  for  each  range  gate.  For  LASBEX,  the  lidar  was  operated  in  the  range-height 
indicator  (RHI)  and  plan  postion  indicator  (PPI)  modes.  In  RHI.  the  beam  of  the  lidar 
sweeps  in  a  ISO0  arc  from  foward  to  overhead  and  then  behind  the  instrument.  In  PPI. 
the  beam  sweeps  at  a  fixed  angle  of  elevation  from  the  horizon  in  a  360°  circle  around 
the  lidar.  the  path  describing  an  inverted  cone.  In  both  modes,  the  display  of  informa- 
tion is  of  the  radial  \  clock)  relative  to  the  lidar.  Familiarity  with  a  doppler  presentation 
is  required  to  analyze  the  PPI  display.  This  analysis  is  similar  to  a  doppler  radar  PPI 
and  is  described  in  Wood  and  Brown  (1983).  In  RHI,  the  lidar  beam  was  aligned  parallel 
to  the  axis  of  the  Salinas  Valley  with  an  orientation  angle  of  140°,  or  east-west  with  an 
angle  of  90°,  with  the  preferred  orientation  being  east-west.  This  was  because 
climatological  records  indicated  that  the  sea  breeze  blew  from  due  west  in  the  Moss 
Landing  area.  In  PPI,  the  elevation  of  the  lidar  beam  was  1°  and  5°.  The  1°  beam 
shows  the  low  level  structure  below  200  meters,  whereas  the  5°  beam  gives  structure 
above  300  meters.  These  values  were  determined  triginometrically.  Table  1  lists  some 
of  the  clear-air  radial  velocity  parameters  of  this  instrument: 


Table   I.     WPL  DOPPLER  RADAR  CAPABILITIES 

Expected  range 

10-15  km 

Range  resolution 

1-40-400  m 

Beam  width 

<  0.1  degree 

Pointing  accuracy 

0.1  degree 

Accuracy  of  radial  velocity  measurment          >   1.0  m  s 

Operation  modes                                          |    PPI,  VAD,  RHI, 

range-time,  raster-scan 

The  lidar  uses  the  backscattered  energy  of  a  laser  from  gas  molecules  and 
particulates  to  detect  atmospheric  features  such  as  clouds  or  aerosol  layers.  Collis  and 
Russell  (1976)  cited  an  example  of  how  lidar  was  used  to  find  the  depth  of  the  mixed 
layer  or  determine  the  heights  of  inversions,  which  can  affect  the  dispersion  of 
pollutants.  Noonkester  et  al.  (1974)  was  able  to  detect  Kelvin-Helmholtz  instabilities  in 
the  boundary-layer.  In  clear  air.  a  lidar  can  detect  the  stratification  associated  with  a 
deep  stable  boundary  layer. 

D.  SURFACE  MEAUREMENTS 

A  portable  coastal  meteorological  station  was  used  to  take  20-second  averages  of 
wind  direction,  wind  speed,  air  temperature,  relative  humidity  and  barometric  pressure 
at  the  surface.  This  instrument  used  a  propeller  anemometer  to  measure  wind  speed. 
A  self-aspirated  thermometer  and  hygrometer  measured  air  temperature  and  relative 
humidity.  Wind  direction  was  taken  from  a  vane  that  gave  angles  with  respect  to  true 
north.  The  ranges  of  the  instrument  errors  were:  ±  5°  for  wind  direction:  +  0.3  m  s  for 
wind  speed;  +  0.3°  for  air  temperature:  and  ±  5%  for  relative  humidity. 

E.  SYNOPTIC  SITUATION  BETWEEN  16-30  SEPTEMBER  1987 

On  16  September.  California  was  under  the  influence  of  a  surface  high  pressure 
system  oil  the  coast  and  a  local  thermal  trough  in  the  center  of  the  state.  The  winds 
were  northerly  to  northwesterly  between  5  and  10  knots.  This  configuration  continued 
to  influence  the  west  coast  from  17  to  24  September.  Fig.  3shows  the  synoptic  situation 
during  this  period.  Around  the  25th,  a  short  wave  which  was  detectable  at  the  850  mb 
level  began  to  pass  through  northern  California  from  the  northwest.  This  wave  built  up 
a  strong  thermal  gradient  off  the  coast  of  California,  providing  an  increase  of  cold  ther- 
mal advection.    This  was  the  most  significant  feature  associated  with  an  actual  synoptic 


frontal  passage  during  LASBEX.  Fig.  4  shows  the  short  wave  location.  Winds  during 
this  wave  passage  were  mainly  northerly  with  speeds  averaging  7  m  s  or  more.  After  the 
passage  of  the  wave  on  the  26th.  high  pressure  off  the  coast  and  the  thermal  trough  in 
the  center  of  the  state  were  reestablished.  This  weather  pattern  persisted  for  the  re- 
mainder of  the   experiment. 

The  combination  o  a  high  pressure  system  over  the  eastern  North  Pacific  Ocean 
and  a  thermal  trough  in  central  California  is  the  typical  summmer  and  early  autumn 
synoptic  pattern.  This  situation  accounts  for  the  common  Monterey  Bay  pattern  of  low 
clouds  and  fog  with  light  and  variable  winds  during  the  night  and  morning.  A  tropical 
disturbance  named  'Otis'  formed  west  of  the  Baja  Peninsula  during  the  23rd  but  did  not 
significantly  affect  surface  and  lower  level  flow  during  that  part  of  the  experiment. 
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Fig.  3.       Typical  Surface  Synoptic  Situation  17-30  September  1987:     Solid  lines  are 

sea  level  isobars.  Values  are  +  1000  mb.  > 


850 


25  September 


Fig.  4.  850mb  Synoptic  Situation  25  September  1987:  Solid  lines  are  isoheights. 
Values  are  +  1000  meters.  Dashed  lines  are  isotherms.  Temperatures  are 
in  decrees  Celsius. 
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III.     DATA  AND  ANALYSIS 

A.     SODAR  THEORY 

1.     Basic  Sodar  Theory  and  Turbulence  in  the  Inertial  Subrange 

// 
This  discussion  of  sodar  theory  and  operation  follows  Lnderwood  (1981).    The 

sodar  utilizes  the  backscatter  of  acoustic  energy  to  determine  certain  parameters  in  the 

atmosphere.    The  relationship  between  the  acoustic  power  emitted  and  the  backscatter- 

ing  cross  section  is  defined  by  the  sodar  equation. 

P(R)  =  P0-4£-Lve~2*R{     .     °  .    V(*)»  P.l) 

'     R~  \   Ao  +  AT  J 

where 

P(R)  =  received  acoustic  power  (watts), 
P0  =  transmitted  acoustic  power  (watts), 
A  =  cross  sectional  area  of  the  receiver  antenna  (tn  ), 
G  =  antenna  gain  factor. 
Ly=  length  of  the  acoustic  pulse  in  space  (m), 
a  =  path  averaged  molecular  attenuation  coefficient  (m~  ), 
R  =  distance  from  the  scatterer  to  the  receiver  (m), 

A0  ..  . 

— —  =     excess    attenuation  and 

A0t  AT 

a(n)  =  backscattering  cross  section  (m~  ). 

The  backscattering  cross  section  can  be  related  to  the  temperature  structure 
parameter  Ot,  and  the  velocity  structure  parameter  Q,  as  shown  by  Little  (1969).  These 
structure  parameters  indicate  the  intensity  of  high  frequency  variations  of  temperature 
or  velocity.    Explicitly,  the  backscatter  cross-section  is 
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c.      2  d  c  - 


(sm4r11/3,  (3.2, 


where 

o{tz)  =  backscattered  acoustic  energy, 

c  =  speed  of  sound, 

T=  absolute  air  temperature, 

Cv  =  velocity  structure  parameter  and 

Cf  =  temperature  structure  parameter. 

Turbulent  mixing  in  clear  air  causes  fluctuations  in  potential  temperature  and 
water  vapor  which  scatter  and  attenuate  the  acoustic  energy  in  the  sodar  signal.  The 
backscattering  cross  section  determines  how  much  energy  is  returned  to  the  sodar  an- 
tenna after  attenuation  by  the  turbulence.  Turbulent  eddies  with  a  size  equal  to  one  half 
the  wavelength  of  the  emitted  acoustic  energy  /.J/2.  (Russell  et  al.,  197-1)  are  detectable 
with  a  sodar.  Eddies  of  this  size  fall  within  the  inertial  subrange  of  turbulence  and  are 
assumed  to  be  isotropic.  Energy  relationships  in  the  inertial  subrange  follow  the 
Kolmogoroff  description  of  isotropic  turbulence.  The  backscattered  data  are  processed 
using  East  Eourier  Transforms  (FFT)  to  establish  the  frequency  shift  due  to  the  radial 
components.  The  two  main  fields  winch  the  monostatic  doppler  sodar  can  map  are  the 
thermal  structure  and  the  mean  velocity  field  along  with  its  variances. 
a.      Thermal  Structure 

Molecular  attenuation  is  described  by  Sutherland  (1975)  and  Piercy  et  al. 
(1977)  as  those  effects  at  the  molecular  scale  which  influence  the  propagation  of  energy 
through  a  medium  such  as  air.  The  factors  which  cause  molecular  attenuation  of  sound 
energy  are  due  to  thermal,  viscous  and  diffusive  effects  as  well  as  rotational  and 
vibrational  losses  on  the  molecular  scale.  The  molecular  attenuation  increases  with  fre- 
quency and  is  the  limiting  factor  in  the  maximum  usable  frequency  in  a  sodar.  The  lower 
limit  is  based  on  antenr  t  size.  "Excess"  attenuation  is  due  to  turbulent  spreading  of  the 
beam  pattern  (Brown  and  Clifford.  1976).  Any  of  the  sodar  signal  that  is  not  attenuated 
can  be  scattered  back  to  rhe  receiver  by  turbulent  eddies  according  to  (3.2). 
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As  the  lapse  rate  changes  in  the  atmosphere,  any  mechanical  turbulence  will 
create  small  eddies  that  backscatter  the  acoustic  energy.  The  monostatic  sodar  (0  = 
180°)  gives  an  indication  of  the  changes  in  the  lapse  rate  along  the  beam  path,  which  can 
be  related  to  thermal  structure  (Wyckoff  et  al.,  1973).  In  a  monostatic  sodar,  only 
temperature  fluctuations  contribute  to  received  power  according  to  (3.2). 
b.      Velocity  Structure 

Doppler  processing  is  used  to  derive  the  velocity  structure.  This  discussion 
follows  Peterson  (19SS)  in  a  handbook  on  doppler  wind  profilers. 

(1)  Basic  Setup  of  A  Doppler  Sodar.  The  sodar  emits  sound  energy  at  a 
narrow  band  of  frequencies  and  is  reflected  by  turbulent  eddies  and  temperature  fluctu- 
ations. A  frequency  shift  (doppler  shift)  is  caused  by  any  motion  on  the  same  axis  as 
the  sodar  beam  and  is  proportional  to  the  speed  of  that  motion.  By  measuring  the 
doppler  shift,  the  radial  velocity  of  the  scattering  volume  can  be  calculated.  In  order  to 
obtain  a  wind  velocity,  Doppler  shifts  from  three  axes  are  required.  The  configuration 
of  the  sodar  is  to  have  one  beam  directed  overhead,  with  the  other  two  at  perpendicular 
cardinal  headings  (i.e.  north  and  east).  The  beams  are  typically  pointed  at  high  ele- 
vations, so  that  they  are  sampling  nearly  the  same  How.  The  vertical  velocity  from  the 
beam  directed  upwards  can  be  subtracted  from  the  other  two.  leaving  the  horizontal 
components  of  the  wind  vector. 

The  sodar  used  during  LASBEX  consists  of  25  separate  transducers. 
Using  electronic  phasing,  these  transducers  were  fired  at  different  times  in  order  to  direct 
the  acoustic  beam.  If  there  was  no  phasing  between  all  the  transducers,  the  beam  would 
travel  directly  o\erhead.  By  delaying  the  time  that  each  transducer  emitted,  the  beam 
of  acoustic  energy  can  be  directed  at  different  angles. 

f2>  Range  Aliasing.  The  signal  from  the  sodar  must  be  of  a  pulse  repe- 
tition frequency  (PRF)  that  precludes  aliasing.  Aliasing  occurs  when  the  the  signal  from 
the  current  pulse  is  overlapped  by  return  signals  from  the  previous  pulse.  These  over- 
lapping signals  will  be  erroneously  interpreted.  This  type  of  aliasing  is  called  range  ali- 
asing because  a  return  from  a  distant  scatterer  could  be  mistaken  for  a  closer  one.  The 
time  between  pulses  must  be  long  enough  to  prevent  this  effect.  The  effective  height  that 
the  sodar  can  function  without  range  aliasing  is  determined  by  the  following  equation: 

Rmax  =2x  PRF  '  (3'3) 

where 
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c  =  speed  of  sound  in  air  (mjs)and 

PRF=  pulse  repetition  frequency  (s-1). 

The  maximum  height  that  the  PSU  sodar  was  set  to  detect  is  750  meters,  which  is  about 
half  its  maximum  range,  and  only  information  starting  50  meters  above  the  surface  was 
recorded.   These  settings  effectively  avoided  any  range  aliasing  during  the  experiment. 

(3)  Velocity  Aliasing.  Only  doppler  shifts  of  certain  frequencies  can  be 
determined  without  error.  Doppler  shifts  that  are  either  too  large  or  too  small  will  be 
incorrectly  mistaken  for  other  velocities.  This  is  called  velocity  aliasing.  The  frequency 
linuts  for  which  a  doppler  shift  can  be  interpreted  correctly  rely  upon  the  sampling  rate. 
The  optimum  selection  is  to  have  a  sampling  rate  as  high  as  possible.  However,  there 
must  be  balance  between  velocity  aliasing  which  requires  a  high  PRF  and  range  aliasing 
winch  requires  a  small  PRF.  The  wind  speed  data  during  LASBEX  showed  no  extreme 
values  and  therefore  velocity  aliasing  was  negligible. 

f4)  Data  Processing.  For  each  detection  cycle,  many  samples  are  gath- 
ered and  converted  into  a  frequency  spectrum  using  the  Fast  Fourier  Transform  tech- 
nique. Doppler  shifts  in  a  frequency  spectrum  show  up  as  peaks  offset  from  0  or  the 
no  shift  position.  Motion  towards  the  sodar  is  indicated  by  a  positive  doppler  shift.  A 
negative  shift  indicates  motion  away  from  the  receiver.  The  doppler  shift  is  converted 
into  a  radial  velocity  component  of  the  wind.  The  height  of  the  frequency  peak  is  a 
measure  of  received  power  and  the  width  of  the  peak  a  measure  of  spread  of  the  acoustic 
energy  over  frequency.  The  spread  of  energy  is  directly  related  to  the  standard  deviation 
of  radial  velocity.  Using  deviations  of  radial  velocities  among  the  three  beams,  the 
standard  deviations  of  each  velocity  component  (au,  al}  oj  can  be  calculated. 

The  preferred  method  is  to  sample  in  the  time  domain,  so  that  several 
consecutive  samples  can  be  gathered  and  averaged  before  conversion  into  the  frequency 
domain  with  an  FFT.  A  mean  frequency  spectrum  is  calculated  from  several  consecutive 
frequency  domains  to  generate  a  spectral  average,  which  gives  a  better  indication  of  a 
doppler  shift  in  the  signal.  The  background  noise  is  filtered  out  to  better  refine  the  fre- 
quency peak.  This  cycle  of  sampling,  conversion,  averaging  and  filtering  is  set  by  the 
user  and  typically  lasts  one  or  two  minutes  per  beam.  The  PSU  sodar  integration  oc- 
curred over  60  seconds,  with  each  beam  sampled  every  10  seconds. 
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B.     LIDAR  THEORY 

The  lidar  uses  the  same  principal  as  the  soJar  to  analyze  the  backscatter  from  a 
pulse  of  energy.  The  lidar  utilizes  a  laser  beam  to  emit  pulses  of  high  intensity  light  that 
will  be  reflected  from  gases,  aerosols  and  particulates.  A  general  form  of  the  lidar 
equation  is  given  (Russell  ct  al.,1974): 

Pr(R)  =  PI^-^Tf](R),2(R),  (3.4) 

""     R 

where 

Pr(R)  =  received  power  from  range,  (R) 

Pt  =  transmitted  optical  power, 

c  =  speed  of  light  in  air, 

r  =  duration  of  laser  pulse. 

A  =  effective  telescope  aperature, 

(J(R)  =  Optical  volume  backscattering  coefficient  and 

t[R)  =  Optical  transmittance  of  atmosphere. 

This  equation  is  similar  in  form  to  the  sodar  equation  (3.1),  except  that  the  lidar 
equation  is  dealing  with  electromagnetic  energy,  whereas  the  sodar  equation  deals  with 
acoustic  energy.  The  received  signal  from  the  lidar  is  a  function  of  the  optical  properties 
and  concentration  of  the  backscattering  volume.  Since  gases  and  aerosols  have  different 
effects  on  the  attenuation  of  light,  the  received  energy  can  be  analyzed  to  determine  what 
causes  the  scatter. 

Aerosols  and  gases  are  distributed  by  the  winds  and  are  trapped  in  layers  due  to  the 
stratification  of  the  atmosphere.  The  detection  of  these  gases  and  aerosols  gives  an  in- 
dication of  the  mean  wind  field  and  the  degree  of  stratification.  The  doppler  lidar  uses 
the  same  principle  of  operation  to  deduce  radial  velocities  as  the  doppler  sodar. 
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C.  SODAR  CROSS  SECTIONS 

1.     Time-Height  Cross  Sections 

The  wind  vector  data  at  the  sodar  were  collected  and  stored  on  floppy  disk  files 
in  tabular  form,  with  a  header  indicating  the  date  and  time  of  the  record.  The  data  were 
recorded  every  25  meters  in  the  vertical,  beginning  50  meters  from  the  surface  up  to  and 
including  750  meters.  These  files  were  transferred  to  the  mainframe  computer  at  the 
Naval  Postgraduate  School  on  mass-storage  for  consolidation  into  separate  files  of  wind 
speed,  wind  direction,  '  ertical  speed,  and  cu,  av  and  aw.  These  consolidated  files  were 
plotted  in  time-height  cross  sections.  Time  for  the  records  was  in  Universal  Time  Co- 
ordinates (UTC),  but  cor  verted  to  decimal  days  for  plotting.  The  conversion  from  UTC 
to  Pacific  Standard  Time  (PST)  is  to  subtract  8  hours  from  UTC.  Plots  were  generated 
in  the  following  formats:  24  hours  up  to  750  meters,  four  hours  up  to  350  meters  and  two 
hours  up  to  550  meters.  The  two  and  four  hour  plots  were  centered  as  close  to  possible 
to  the  passage  of  the  sea  breeze  front.  The  24-hour  plots  started  and  ended  at  0600  local 
time  (1400  UTC  or  0.5833  decimal  day).  The  raw  data  from  mass-storage  were  read  and 
set  to  a  standard  grid  with  a  resolution  of  five  minutes  on  the  abscissa  and  25  meters 
on  the  ordinate.  The  density  of  the  grid  depended  on  the  time  and  height  resolution  of 
the  desired  field  (i.e.  24-hour  plots  with  750  meter  heights  were  29  by  288,  whereas  two 
hour  plots  with  550  meter  heights  were  19  by  24).  Contours  were  fitted  to  the  standard 
grid  using  a  rational  spline  interpolation  scheme.  Plots  were  generated  and  checked  for 
bad  data,  which  showed  up  as  unusually  strong  gradients  surrounding  a  single  point. 
Vertical  records  were  manually  checked  for  consistency  and  errant  readings  were  re- 
moved by  comparison  to  adjacent  values.  The  two  hour  cross-sections  of  wind  direction, 
wind  speed,  vertical  velocity  and  g„  are  included  as  appendices  A,  B.  C  and  D.  The 
24-hour  cross-sections  of  wind  direction  and  wind  speed  are  included  as  appendices  U 
and  F. 

D.  DERIVED  DATA 

Horizontal  divergence  and  turbulent  kinetic  energy  were  derived  from  the  velocity 
and  o,_.  o,  and  a„  fields  respectively  and  plotted  in  the  same  format  as  the  velocity  fields. 

Horizontal  divergence  was  determined  from  the  continuity  equation  for 
incompressible  flow: 

ilL  +  4L  +  ZlL\=0.  (3.5) 

ex        ov        c: 


is 


The  divergence  was  calculated  from  the  vertical  velocity  field  such  that 

4^  —  v^.r.  (3.6) 

This  calculation  was  simple,  since  the  vertical  velocity  component  of  the  wind   with   re- 
spect to  height  was  readily  available  from  the  sodar  data. 

The  turbulent  kinetic  was  determined  using  the  values  of  ctu,  ou  and  aw  as  follows: 

e  =  -i-  (<rj  +  a\  4-  ol),  (3.7) 

where  e  equals  turbulent  kinetic  energy. 

E.     ANALYSIS  OF  TIME-HEIGHT  CROSS  SECTIONS 
1.     Wind  Direction 

a.     Wind  Direction  Gradients  and  Onset  of  Sea  Breeze 

The  criteria  that  mark  the  onset  of  the  sea  breeze  typically  are  the  change 
in  wind  direction  and  increase  of  wind  speed.  The  wind  data  from  the  sodar  were  ana- 
lyzed to  determine  the  degree  of  the  wind  speed  change.  A  sharp  change  in  wind  direc- 
tion with  time  which  indicated  sea  breeze  frontal  passage  was  present  in  10  plots  during 
LASBEX.  Wind  direction  was  contoured  with  an  interval  of  30°.  Onshore  wind  was 
considered  270°  ±  60°  and  offshore  90°  ±  60°.  Up  valley  was  330  ±  20°.  Down  valley 
was  considered  150  +  20°.  Frontal  passage  and  the  onset  of  the  sea  breeze  was  evident 
when  the  wind  direction  showed  a  strong  temporal  gradient.  The  mean  temporal  gradi- 
ent for  sea  breezes  during  LASBEX  was  on  the  order  of  a  90°  clockwise  shift  over  44 
minutes  for  a  strong  offshore  ilow  and  o\er  42  minutes  for  a  weak  offshore  flow.  The 
fastest  change  occurred  over  10  minutes  during  strong  offshore  flow,  while  the  slowest 
change  occurred  over  104  minutes  during  weak  offshore  flow.  90°  was  selected  as  the 
criterion  to  define  sea  breeze  onset,  when  the  wind  turned  through  alongshore  to 
onshore  flow,  because  the  observed  temporal  wind  shift  was  greatest  over  this  interval. 
The  relative  strength  of  the  offshore  flows  was  divided  into  strong  and  weak  flows  based 
on  the  preexisting  winds  below  400  meters  before  the  onset  of  the  sea  breeze.  The 
400-meter  level  was  selected  since  it  was  the  average  depth  of  the  established  offshore 
flow.  Strong  flows  were  considered  3  m/s  or  greater  and  weak  flows  less  than  3  m  s. 
There  was  no  correlation  between  the  sharpness  of  the  front  based  on  the  temporal 
gradient  and  the  strength  of  the  previous  oifshore  flow.    This    agrees    with    the  model 
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result  of  Pearson  et  al.  (1983)  that  the  intensity  of  the  sea  breeze  circulation  is  unaltered 
by  the  mean  flow.  Table  2  on  page  20  lists  the  various  gradients  for  sea  breeze  fronts 
that  were  detected  during  LASBEX. 


Table  2.  GRADIENTS  OF  WIND  DIRECTION  AT  THE  SEA  BREEZE 
FRONT:  \is  =  mean  of  strong  offshore  flows,  /u„  =  mean  of  weak  off- 
shore flows 


Date 

Temporal  Direc- 
tional V(°  min) 

Vert.  Directional 
V(°/m) 

Oft  shore  Flow  (Strong 
or  Weak) 

16 

90  24 

180  45 

Strong 

17 

90  14 

ISO  25 

Weak 

IS 

90  70 

ISO  37 

Weak 

19 

90  10 

ISO  50 

Strong 

20 

90  60 

ISO  105 

Strong 

21 

90  34 

180/45 

Strong 

22 

90  30 

1  SO  90 

Weak 

23 

90  104 

ISO  ISO 

Weak 

77 

90  82 

ISO  50 

Strong 

29 

90  56 

ISO  50 

Strong 

f-3 

90  44 

ISO  58 

n  a 

(J.w 

90  42 

ISO  S3 

n  a 

A  shear  layer  was  evident  in  the  direction  plots  as  a  strong  gradient  in  the 
vertical.  This  shear  zone  indicates  the  boundary  between  the  offshore  flow  and  the  sea 
breeze.  These  gradients  are  shown  in  Table  2.  The  mean  wind  shift  in  the  vertical  was 
18o°  over  58  meters  and  S3  meters  for  the  strong  and  weak  offshore  cases,  respectively. 
It  does  appear  that  a  strong  offshore  flow  limited  the  vertical  extent  of  the  shear  zone 
and  that  a  weak  flow  allowed  a  deeper  region  of  directional  shear,  which  agrees  with  the 
results  produced  by  Simpson  and  Britter  (1980)  in  a  laboratory  experiment,  as  shown  in 
Fig.  5  on  page  21. 

In  the  majority  of  the  frontal  passages,  the  wind  direction  change  preceeded 
a  change  in  wind  speed  or  vertical  velocity.  In  all  cases,  the  direction  of  the  breeze  after 
the  direction  shift  was  300°  ±  20°, which  was  onshore  and  up  valley. 

The  average  onset  times  for  the  sea  breeze  determined  at  the  sodar  site  by 
the  change  in  wind  direction  was  1742  UTC.    The  earliest  onset  time  was  at  1644  L'TC 
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Fig.  5.        Photographs    of   a    Laboratory    Gravity    Head    Structure    in    Prevailing 

Flows:       The  top  photo  is  taken  during  opposing  flow.  The  bottom  is 
taken  during  accompanying  flow,  (from  Simpson  and  Britter,  1980) 

and  the  latest  was  at  ISjO  UTC.   The  time  for  the  increase  in  wind  speed  and  in  vertical 
velocity  were  also  tabulated  from  the  sodar  data. 

Onset  times  were  also  determined  from  the  surface  data.  Onset  was  indi- 
cated by  a  change  in  wind  direction  and  speed,  temperature  and  relative  humidity. 
Shortly  before  sea  breeze  onset,  the  wind  direction  was  very  variable  with  the  trend 
changing  from  ollshore  to  onshore  flow.    Wind  speeds  before  onset  decreased  to  near 
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zero,  then  rapidly  increased  as  the  sea  breeze  front  passed.  Air  temperature  gradually 
rose  due  to  solar  heating  in  the  morning,  peaked  at  frontal  passage  and  then  decreased 
as  the  sea  breeze  was  established.  Relative  humidity  at  onset  increased  quickly  with 
frontal  passage  and  leveled  off  at  the  sea  breeze  continued.  A  typical  surface  record  with 
sea  breeze  passage  is  shown  in  Fig.  6. 

In  most  cases,  the  wind  direction  shift  was  consistent  with  the  observed 
sodar  wind  change.  The  mean  time  of  onset  measured  by  the  surface  data  was  1813 
UTC.  The  earliest  frontal  passage  occurred  at  1715  UTC.  while  the  latest  passage  oc- 
curred at  1945  UTC.  The  variability  of  wind  direction  shifts  in  the  surface  measurements 
and  sodar  as  well  as  the  absence  of  data  on  various  days  from  both  sources  explains  the 
difference  between  the  surface  and  sodar  mean  values.  In  Table  3,  the  times  of  onset 
based  on  the  surface  data  at  the  sodar  site  and  the  directional  shift  are  shown,  as  well 
as  the  times  when  the  speed  and  vertical  velocity  increases  occurred. 


Table  3.  TIMES  OF  ONSET  FOR  THE  SEA  BREEZE:  sfcso  =  Onset  based  on 
Surface  Measurements  at  the  Sodar.  sdir=  Onset  of  wind  direction  shift 
at  Sodar.  sspd=  speed  change  at  Sodar,  svrt=Time  of  Maximum  Updrafts 
at  Sodar.  All  times  in  UTC. 


Date 

sfeso 

sdir 

sspd 

svrt 

16 

n  a 

1745 

1750 

1750 

17 

1715 

1714 

1717 

1715 

IS 

1 750 

1748 

1813 

1 750 

19 

1730 

1 8oo 

1819 

1820 

2<) 

n  a 

1644 

1709 

1659 

21 

n  a 

1729 

1749 

1  ~5o 

22 

n  a 

1715 

1730 

1730 

23 

1830 

1830 

1850 

1850 

->s 

1945 

n  a 

n  a 

n  a 

26 

1S55 

n  a 

n'a 

n  a 

27 

1830 

1818 

1824 

1819 

29 

1800 

1  756 

1801 

1801 

3() 

1745 

n  a 

n  a 

n  a 

M 

1813 

1742 

1757 

1752 

22 
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Fig.  6.      Surface  Data  Record  from  29  September  1987  at  the  Sodar  Site 

b.     Height  and  Duration  of  the  Sea  Breeze 

The  directional  gradient  which  indicates  the  boundary  between  the  onshore 
and  offshore  flow  could  be  followed  in  the  24-hour  contour  plots  as  the  sea  breeze 
deepened.  The  maximum  vertical  extent  of  the  sea  breeze  was  determined  for  the  10 
cases  and  had  a  mean  elevation  of  659  meters,  and  ranged  between  570  and  800  meters.' 

The  cessation  time  of  the  sea  breeze  was  indicated  by  a  counter-clockwise 
turning  of  the  winds,  a  decrease  in  the  height  of  the  sea  breeze  and  the  establishment 
of  offshore  flow.  Duration  of  the  sea  breeze  was  the  difference  between  the  onset  and 
cessation  times  determined  from  the  directional  data.  The  mean  duration  of  the  sea 
breeze  was  10.7  hours.  The  longest  duration  was  13.4  hours  and  the  shortest  was  8.1 
hours.  The  mean  duration  for  the  sea  breeze  was  also  determined  using  surface  data 
measurements  at  the  lidar  and  sodar  sites.    Using  the  surface  data,  cessation  was  when 
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the  wind  direction  was  alongshore  and  the  wind  speed  had  dropped  to  2  in  s.  The  mean 
durations  calculated  using  surface  data  at  the  sodar  and  lidar  were  10.2  and  11.4  hours 
respectively.  A  summary  of  maximum  height  and  duration  of  the  sea  breeze  is  given  in 
Table  4. 


Table 

4.     MAXIMUM  HEIGHT  AND  DURATION  OF  THE  SEA  BREEZE 

Date 

Maximum 
Height  (meters) 

Duration  by 
Sodar  (h) 

Duration  by 

Surface  at  Sodar 

(h) 

Duration  by 

Surface  at  Lidar 

(h) 

16 

600 

S.l 

n  a 

n  a 

17 

725 

9 

13.0 

n'a 

IS 

640 

13.4 

13.0 

12.0 

19 

SOU  (est) 

13.0 

n  a 

n  a 

20 

600 

11.3 

n  a 

11 

21 

570 

n'a 

n'a 

10 

•)■} 

680 

11.5 

n  a 

n  a 

7-5 

650 

8.5 

12.0 

12.0 

■>^ 

n  a 

n  a 

10.0 

12.0 

26 

n  a 

n  a 

n  a 

8 

27 

625 

11. S 

11 

9 

28 

n  a 

n  a 

n  a 

8 

29 

7()0 

12.2 

10 

10 

pL 

659 

10.7 

11.4 

10.2 

In  Fig.  7  on  page  25,  the  increasing  height  of  the  sea  breeze  is  clearly  shown,  as  well  as 
the  maximum  height  achieved  and  the  decrease  in  the  height  with.  time. 
c.     Slope  of  Sea  Breeze  Fronts 

The  slope  of  the  sea  breeze  front  was  determined  for  the  first  kilometer 
horizontally  of  the  front  and  for  the  established  sea  breeze  flow  subsequent  to  the  first 
kilometer.  The  slope  for  the  first  kilometer  of  the  front  was  calculated  by  measuring  the 
difference  from  the  50  meter  level  to  the  highest  height  of  the  onshore  flow  in  the  first 
kilometer  and  dividing  by  the  horizontal  distance  between  these  points.  The  slope  for 
the  rest  of  the  sea  breeze  was  calculated  by  measuring  the  difference  between  the  greatest 
height  in  the  first  kilometer  and  the  greatest  height  of  the  onshore  flow  during  rest  of  the 
sea  breeze.     Horizontal  distance  was  calculated  by  multiplying  the  time  difference  be- 
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LASBEX  WIND  DIRECTION  CROSS  SECTION 
SEPTEMBER  1987 


m 


m  h 


wMi^A  tot 

m  m  ifl  m  y 


16.708  16.833  16.958  17.083  17.208  17.333  17.458  17.583 

TIME(DECIMAL  DAY) 

CONTOUR  INTERVAL:  30  DEGREES 


Fig.  7.      24  Hour  Directional  Plot  for  16  September  1987 

tween  the  two  heights  by  the  speed  of  advance  of  the  front.    Speeds  of  advance  were 
found  by  taking  the  time  difference  between  time  of  frontal  passage  based  on  the  surface 
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data  at  the  lidar  and  the  time  of  passage  at  the  sodar.  For  those  cases  where  a  speed 
of  advance  could  hot  be  calculated,  the  mean  speed  of  advance  was  used.  The  mean 
speed  of  advance  during  LASBEX  was  2.5  m/s.   These  speeds  are  listed  in  Table  5. 


Table  5.  SPEED  OF  ADVANCE  OF  THE  SEA  BREEZE  FRONT:  Times  of 
frontal  passage  are  from  the  surface  data,  except  in  those  marked  *,  which 
are  based  on  sodar  data. 


Date 

Lidar  (ETC) 

Sodar  (ETC) 

At(h) 

Speed  of  Ad- 
vance (m,'s) 

IS 

1 720 

1750 

0.5 

3.1 

19 

1 700 

1730 

0.5 

3.1 

20 

1530 

1644* 

1.5 

1.0 

21 

1 700 

1729- 

.48 

3.2 

23 

1 700 

IS  30 

1.5 

1.0 

25 

1915 

19-45 

0.5 

3.1 

26 

1815 

1S55 

0.67 

2.3 

27 

1 800 

IS  30 

0.5 

3.1 

29 

1708 

1750 

0.7 

-i  -> 

30 

1730 

1  soo 

0.5 

3.1 

u 

1724 

1806 

0.74 

2.5 

In  the  first  kilometer,  the  mean  slope  was  1:12  for  strong  offshore  flow  and 
1:10  for  weak  flow.  The  typical  slope  after  the  first  kilometer  was  1:94  for  strong  off- 
shore flow  and  1:160  for  weak  offshore  flow.  All  the  values  except  the  subsequent  slope 
for  the  weak  cases  are  similiar  to  ones  cited  in  Simpson  et  al.  (1977),  which  were  1:4  and 
1:10  for  the  first  kilometer  and  1:100  for  the  subsequent  portion  of  the  sea  breeze. 
However,  there  were  no  slopes  on  the  order  of  the  one  reported  by  Helmis  et  al.  (19S7;, 
which  had  a  slope  of  1:3  in  the  first  kilometer.  There  appears  to  be  no  correlation  be- 
tween the  degree  of  slope  in  the  first  kilometer  and  subsequent  to  passage.  That  is,  a 
strong  slope  in  the  first  kilometer  did  not  mean  that  the  rest  of  the  sea  breeze  also  had 
a  sharp  slope,  and  vice  versa.  This  indicated  that  the  sea  breeze  was  very  dynamic  and 
its  velocity  and  structure  changed  throughout  its  duration.  The  slopes  for  the  sea  breeze 
fronts  during  LASBEX  are  summarized  in  Table  6. 
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Table  6. 


SLOPE  OF  THE  SEA  BREEZE  FRONT: 

flows.  u„  =  mean  of  weak  offshore  flows. 


A*.  =  mean  of  strong  offshore 


Date 

Slope  in  first  kilome- 
ter 

Slope  after  first  kilo- 
meter 

Offshore  Flow 
(Strong  or  Weak) 

16 

1:17 

1:87 

Strong 

17 

1:5 

1:158 

Weak 

IS 

1:8 

1:205 

Weak 

19 

1:6 

1:153 

Strong 

20 

1:11 

1:110 

Strong 

21 

1:9 

1:59 

Strong 

22 

1:14 

1:58 

Weak 

23 

1:11 

1:220 

Weak 

27 

1:15 

1:100 

Strong 

29 

1:16 

1:56 

Strong 

u: 

1:12 

1:94 

n  a 

jU„ 

1:10 

1:160 

n  a 

2.     Wind  Speed 

a.     Synopsis  of  H  'hid  Speed  Data  During  LAS  BEX 

A  characteristic  of  the  sea  breeze  that  was  investigated  during  this  study  is 
the  difference  in  time  between  the  change  of  wind  direction  and  the  increase  in  wind 
speed  and  whether  or  not  this  time  difference  was  affected  by  the  previous  offshore  flow. 
The  general  structure  of  the  wind  speed  field  was  also  compared  with  the  wind  direction 
field  to  see  whether  features,  such  as  shear  zones,  indicated  in  one  agreed  dimensionally 
with  the  other,  that  is  if  the  height  and  width  of  the  shear  zone  in  one  case  was  similar 
to  the  other. 

The  wind  speed  data  were  plotted  with  contour  intervals  of  two  meters  per 
second  for  the  24  hour  plots  and  one  meter  per  second  for  the  shorter  plots.  These  plots 
were  analyzed  for  a  continuous  increase  of  wind  speed  with  time.  The  mean  time  of  wind 
speed  increase  at  the  sodar  site  was  1757  UTC  or  17  minutes  after  the  shift  in  wind  di- 
rection. The  mean  temporal  gradient  was  a  1  m.'s  change  over  55  minutes  for  strong 
offshore  flow  and  1  m's  over  21  minutes  for  a  weak  offshore  flow.  In  this  case,  the 
findings  do  not  agree  with  Pearson's  (1983)  model  that  shows  no  relation  between  the 
intensity  of  the  mean  flow  and  the  intensify  o[  the  sea  breeze  circulation.    The  results 
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from  LASBEX  point  to  the  shortcomings  that  a  model  has  in  trying  to  represent  a  dy- 
namic How  such  as  the  sea  breeze  circulation. 

There  was  a  shear  zone  at  the  boundary  between  the  sea  breeze  and  the 
oil  shore  flow.  The  vertical  gradient  of  this  shear  zone  was  a  1  m  s  change  over  40  meters 
for  the  strong  case  and  1  m/s  over  44  meters  for  the  weak  offshore  flows.  The  level  of 
the  wind  speed  shear  correlated  with  the  level  of  the  wind  direction  shear.  Table  7  lists 
the  gradients  in  speed. 


Table  7. 


GRADIENTS  OF  WIND  SPEED  AT  THE  SEA 

=  mean  of  strong  offshore  flow,  fxs=  mean  of  we 


BREEZE  FRONT:      iui 

ak  offshore  flows. 


Date 

Temporal  Speed 
V(m  s  min) 

Vert.  Speed  V(m,  s  m) 

Offshore  Flow  (Strong 
or  Weak) 

16 

1/96 

1  40 

Strong 

17 

1  22 

1  46 

Weak 

18 

1  40 

1  50 

Weak 

19 

1  34 

1  '35 

Strong 

20 

1  100 

1  24 

Strong 

21 

1  66 

1  36 

Strong 

7-) 

1  15 

1  20 

Weak 

23 

1  6 

1  58 

Weak 

77 

1  21 

1  55 

Strong 

29 

1   10 

1  50 

Strong 

u. 

1  55 

1  40 

n  a 

flu 

F2J 

1  44 

n  a 

3.     Wind  speed  aftei  frontal  passage 

The  mean  wind  speed  within  three  hours  after  frontal  passage  during  LASBEX 
was  5.3  m/s  for  strong  offshore  flows  and  ranged  between  4  and  6  m/s.  The  mean  for 
weak  offshore  flow  was  7  m/s  and  was  between  6  and  10  m  s.  The  wind  speeds  of  the 
sea  breeze  reached  maximum  intensity  about  the  same  time  as  maximum  height  was 
achieved^  The  mean  maximum  speeds  during  the  entire  duration  of  the  sea  were  8.3  m  s 
for  stong  flows  and  7  m/s  for  weak  flows.  The  maximum  speeds  ranged  between  6  and 
10  m/s.  Offshore  flow  was  clearly  present  before  the  onset  of  the  sea  breeze.  This  flow 
was  considered  the  established  land  breeze  and  was  usually  confined  below  400  meters. 
The  mean  speed  of  the  strong  offshore  flow  was    6.3  m  s    and    the    mean  of  the  weak 
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oflshore  flows  was  2.0  m  s.  All  the  weak  offshore  flows  were  2  m  s  or  less,  while  the 
strong  offshore  flows  were  between  4  and  S  m  s.  Table  8  on  page  29  gives  the  various 
maxima  of  wind  speed. 


Table  8.  MAXIMUM  WIND  SPEEDS  DURING  THE  SEA  BREEZE  AND 
MAXIMUM  OFFSHORE  FLOW:  ^  =  mean  of  strong  oflshore  flows, 
f.mu  sub  w  =  mean  of  weak  offshore  flows. 


Date 

Maximum  Speed 
after  onset  (m's) 

Maximum  dur- 
ing Sea  Breeze 
(m  s) 

Maximum  Previ- 
ous Offshore 
Flow  (m  s) 

Olfshore  flow 

(Strons  or 

Weak) 

16 

6 

10 

8 

Strong 

17 

6 

6 

2 

Weak 

IS 

6 

6 

2 

Weak 

19 

4 

10 

6 

Strong 

2(» 

6 

10 

6 

Strong 

21 

4 

6 

4 

Strong 

22 

10 

10 

■> 

Weak 

2 1 

6 

6 

-> 

Weak 

2/ 

6 

6 

8 

Strong 

29 

6 

10 

6 

Strong 

,/.;, 

5.3 

8.7 

6.3 

ii  a 

/U„ 

7 

- 

2 

n  a 

As  the  sea  breeze  established  itself,  the  top  of  the  onshore  breeze  increased  until 
it  assumed  a  relatively  constant  elevation.  Above  the  this  level,  the  return  flow  which 
completes  the  sea  breeze  circulation  was  detectable  by  the  sodar.  Since  the  sodar  signal 
had  little  data  past  750  meters,  the  extent  and  characteristics  of  the  return  flow  could 
not  be  determined.  However,  the  information  from  the  lidar  clearly  showed  oflshore 
flow  on  16  September  that  began  around  400  meters  and  extended  vertically  up  to  1600 
meters. 

4.     Vertical  Velocit  •  Structure 

a.      Vertical  Velocity  Maxima  and  Location  of  the  Updraft  Region 

The  vertical  velocity  was  investigated  in  order  to  determine  the  strength  and 
location  of  the  updraft  region  in  the  sea  breeze  front.  Vertical  velocity  cross-sections 
showed  good  correlation  with  both  the  wind  direction  and  speed  cross-sections.  The 
maximum  in  upward  vertical  velocity  occurred  after  the  wind  direction  change  and  about 
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the  same  time  as  the  increase  in  wind  speed  (see  Table  3  on  page  22).  The  occurrence 
of  the  updraft  region  indicated  the  front  of  the  longitudinal  vortex  which  is  the  major 
current  within  the  sea  breeze  frontal  circulation.  The  mean  maximum  vertical  velocity 
observed  after  frontal  passage  was  0.56  m/s,  with  the  range  of  velocities  between  0.2  and 
1.0  m/s.  This  vertical  velocity  is  much  smaller  than  the  maximum  velocity  observed  by 
Ogawa  et  al.  (1986)  of  3.0  m/s;  but  similar  to  the  values  reported  by  Simpson  et  al. 
(1977)  of  0.7  to  3.0  m/s  and  the  average  value  that  Helmis  et  al.  (1987)  observed  of  0.5 
m/s.  A  summary  of  the  vertical  characteristics  of  the  sea  breeze  is  contained  in 
Table  9. 


Table 

9.     VERTICAL  MOTION  CHARACTERISTICS  OF  THE  SEA  BREEZE 

Date 

Max.  Vert.  Ex- 
tent (m) 

Max.  Vert.  Ve- 
locity (m/s) 

Location  o[ 
Max.  Velocity 

Offshore  Flow 

(Strong  or 

Weak) 

10 

225 

0.6 

50-150 

Strong 

17 

350 

0.4 

100-275 

Weak 

IS 

370 

1.0 

50-250 

Weak 

19 

270 

0.4 

50-200 

Strong 

20 

1 50 

0.6 

50-175 

Strong 

21 

200 

0.6 

50-175 

Strong 

2 1 

3  3  5 

0.6 

75-225 

Weak 

J^ 

3  SO 

0.8 

75-350 

Weak 

27 

220 

0.6 

50- 1 50 

Strong 

29 

175 

0.2 

50-150 

Strong 

u 

2<)7 

0.5 

50-167 

n  a 

ult 

359 

0.7 

75-275 

n  a 

The  location  of  the  vertical  velocity  maximum  within  the  sea  breeze 
front  gravity  head  structure  is  shown  in  Fig.  8,  where  the  longitudinal  \ortex  turns  up- 
ward. 

The  location  of  the  updraft  region  was  found  to  be  related  to  the  height  of 
the  sea  breeze  front,  which  in  turn  was  related  to  the  strength  of  the  previous  offshore 
flow.  The  mean  location  of  the  maximum  vertical  velocity  was  at  207  meters  in  strong 
onshore  flow  and  359  meters  in  weak  offshore  flows.  When  offshore  winds  were  strong 
before  the  sea  breeze,  the  updraft  region  started  at  a  lower  elevation  and  died  out  at  an 
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Fig.  8.  Flow  regions  of  the  Gravity  Head:  The  stippled  region  is  where  the  op- 
posing flow  is  being  overrun  by  the  gravity  head  current  (from  Simpson 
and  Britter,1980) 

average  of  of  1S3  meters,  and  when  the  winds  were  weak,  the  updraft  began  around  75 
meters  and  died  out  around  250  meters.  These  conditions  both  agree  with  the  gravity 
head  structure  in  prevailing  flow  that  was  observed  by  Simpson  and  Britter  (1980).  In 
a  strong  opposing  current,  the  gravity  head  is  flatter  with  the  flow  constrained  to  a 
shallower  region  than  in  onshore  or  weak,  offshore  flow.  The  leading  edge  of  the  onshore 
flow  also  tends  to  have  a  more  pronounced  "knee"  where  it  is  overriding  the  counter- 
How,  (see  Fig.  5  on  page  21)  The  updraft  in  the  onshore  flow  is  higher,  indicating  some 
influence  from  the  counter-flow.  Fig.  S  shows  how  the  updraft  is  affected  by  the  over- 
running and  elevated  region  at  the  front  of  the  gravity  head  flow.  The  "knee"  structure 
of  the  gravity  head  current  was  apparent  in  all  the  sea  breeze  fronts  investigated  during 
LASBEX. 
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5.     Variations  of  the  Vertical  Velocity 

The  standard  deviation  of  the  vertical  velocity  was  used  to  determine  those  areas 
of  the  sea  breeze  front  that  were  the  most  turbulent.  The  standard  deviation  is  higher 
in  turbulent  regions  and  near  zero  in  non-turbulent  stable  layers.  The  values  for  o„  were 
directly  generated  by  the  PSU  sodar.  These  values  were  read  and  plotted  in  the  same 
manner  as  the  other  parameters.  Generally,  the  values  of  a„  decreased  with  height, 
which  indicated  a  reduction  in  the  amount  of  turbulence.  Vertical  velocity  fluctuations 
were  damped  by  the  boundary  at  the  base  of  the  inversion.  Table  10  on  page  33  shows 
the  change  of  ow  across  the  front  at  80  meters. 

The  values  of  ow  did  not  consistently  decrease  with  time  at  the  lower  levels  (i.e 
<  100  m)  after  the  passage  of  the  sea  breeze  front.  In  the  majority  of  the  cases,  the 
value  of  a„  increased  by  an  average  of  0.17  m  s,  with  other  maxima  occurring  at  various 
time  in  the  established  sea  breeze  flow.  The  changes  in  ow  were  between  -0.1  and  +  0.4 
m,  s.  Helmis  et  al.  (1987)  observed  that  the  values  of  aw  decreased  after  the  passage  ol' 
the  sea  breeze,  whereas.  Ogawa  et  al.  (1986).  found  that  the  values  increased  across  the 
front.  The  explanation  given  by  Helmis  was  that  the  times  of  onset  and  the  location  are 
the  governing  factors  in  the  degree  of  turbulence  generated  by  the  sea  breeze.  In  the 
case  of  LASBEX,  the  onset  of  the  sea  breeze  occurred  early  in  the  day.  The  variation 
in  vertical  velocity  was  mainly  due  to  mechanical  mixing  at  and  within  the  front.  Helmis 
observed  sea  breezes  that  passed  in  the  early  afternoon.  The  sea  breeze  in  this  case  acted 
to  decrease  the  amount  of  thermal  activity  caused  by  the  heated  land.  There  appeared 
to  be  no  correlation  between  the  strength  of  the  previous  offshore  flow  and  the  change 
of  <7„  across  the  front.  In  the  case  study  for  16  September,  the  behavior  of  ou  will  be 
investigated. 

F.     DERIVED  DATA  TIME-HEIGHT  CROSS  SECTIONS 

1.     Horizontal  Divergence 

Since  the  sea  breeze  front  is  an  area  of  strong  convergence,  it  was  logical  to 
calculate  the  intensity  of  the  convergence  from  the  velocity  fields.  The  horizontal  di- 
vergence was  calculated  using  the  vertical  velocity  field  according  to  equation  (3.6).  The 
divergence  calculation  using  the  vertical  velocity  field  was  straight  forward,  since  the 
sodar  produced  velocities  as  a  function  of  z. 
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Table  10. 


STANDARD    DEVIATION    OF    VERTICAL    VELOCITY    ACROSS 
THE  SEA  BREEZE  FRONT 


Date 

oH  Before 

a„  After 

Increase  or  De- 
crease 

Oilshore  Flow 

(Stronc  or 

Weak) 

16 

0.9 

o.s 

1 

Strong 

17 

o.s 

0.8 

no  change 

Weak 

18 

0.9 

1.3 

T 

Weak 

19 

1.0 

0.9 

I 

Strong 

2o 

0.7 

0.9 

T 

Strong 

21 

l.o 

0.8 

1 

Strong 

2"> 

0.7 

0.8 

T 

Weak 

23 

0.7 

0.8 

t 

Weak 

27 

0.8 

0.9 

T 

Strong 

29 

0.6 

0.7 

T 

Strong 

R 

0.73 

0.90 

0.17 

n  a 

,u„ 

0.97 

0.83 

0.14 

n  a 

Lyons  and  Oisson  (1973)  calculated  and  plotted  horizontal  divergence  in  time- 
distance  coordinates,  using  the  u-component  of  velocity  based  on  the  assumption  that 
all  variations  occurred  along  a  two-dimensional  axis  (i.e.  dvjdy  =  0),  according  to  the 
continuity  equation.  The  maximum  values  for  convergence  that  they  calculated  at  the 
lake  breeze  front  was  on  the  order  of  2  x  10_35-1.  The  values  of  convergence  calculated 
during  LASBLX  exceeded  A  x  10_3s-1  at  the  head  of  the  front,  and  at  times  were  as  much 
as  12  x  10_35_1.  The  values  calculated  from  the  vertical  velocity  field  of  the  sodar  are 
believed  to  be  a  better  representation  of  the  actual  divergence  along  sea  breeze  front. 
because  the  values  of  vertical  velocity  were  directly  measured.  The  lake  breeze  studied 
by  Lyons  and  Oisson  had  lower  computed  values  of  maximum  vertical  velocity,  around 
0.105  to  0.130  m  s,  whereas  during  LASBEX,  maximum  vertical  velocities  were  almost 
an  order  of  magnitude  larger  (i.e.  O.S  to  1.0  m.'s).  This  is  further  evidence  that  Lyons' 
and  Olsson's  calculations  may  have  been  low  and  that  the  assumption  of  two- 
dimensional  flow  was  perhaps  incorrect.  Fig.  9  on  page  35  shows  the  typical  divergence 
plots  from  Lyons  and  Oisson.  Fig.  10  on  page  36  is  an  example  of  the  horizontal  di- 
vergence calculated  from  the  vertical  velocity  field  on  20  September  19S7. 
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Fig.  9.  Divergence  values  of  the  Chicago  Lake  Breeze  Study:  Contour  interval  is 
50  x  JO'-V.  Shaded  areas  indicate  convergence,  (from  Lyons  and  Olsson, 
1973) 
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LASBEX  HORIZONTAL  DIVERGENCE  CROSS-SECTION 
SEPTEMBER  1987 


20  6654  20  7071  20  7488  20  7905  20 

TIME  (DECIMAL  DAY) 
4  HR  PLOT  OF  SEA  BREEZE  CONTOUR  INTERVAL:  4E-03/SECOND 


Fig.   10.      Divergence  values  of  the  LASBEX  Study:     Convergence  at  the  front  oc- 
curs near  20.7071  decimal  day. 
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2.     Turbulent  Kinetic  Energy  (TKE) 

Since  turbulence  and  mechanical  mixing  is  highest  within  the  head  of  the  sea 
breeze  front,  the  locations  of  maximum  TKE  should  coincide  with  those  regions  of 
strong  shear,  strong  convergence,  or  overturning  flow  that  creates  turbulent  eddies.  A 
general  view  of  these  regions  was  conducted  using  the  standard  deviations  from  the  ve- 
locity fields.  The  contours  for  TKE  were  calculated  by  squaring  the  standard  deviations 
from  the  three  velocity  components  obtain  their  variance  for  substitution  in  equation 
(3.7).  The  values  of  TKE  were  then  plotted  in  two  hour  format.  TKE  values  were  found 
to  be  highest  at  the  base  and  top  of  the  gravity  head,  where  the  longitudinal  vortex  is 
overturning  and  interacting  with  both  the  opposing  flow  and  the  vertical  shear  region 
above  the  gravity  head,  corresponding  to  the  area  where  horizontal  flow  is  directed  up- 
wards and  into  the  shear  zone  within  the  gravity  head.  The  case  study  for  16  September 
will  show  these  characteristics. 
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IV.     CASE  STUDY  OF  THE  SEA  BREEZE  ON  16  SEPTEMBER  1987 

A.  MOTIVATION  FOR  STUDY 

'September  16,  19S7  proved  to  be  an  unusual  day  since  the  sodar  detected  the  pas- 
sage of  two  sea  breeze  fronts.  The  sea  breeze  front  that  passed  in  the  morning  is  the 
subject  of  this  case  study.  The  late  sea  breeze  is  currently  under  further  investigation. 
Since  the  16th  was  unusual,  WPL  provided  the  lidar  information  gathered  during  the  two 
sea  breeze  passages.  Surface  observations  and  radiosonde  observations  were  not  avail- 
able, but  were  deemed  unnecessary  in  the  comparison  of  the  wind  field  structure  of  the 
sea  breeze  that  was  detected  by  both  the  lidar  and  the  sodar. 

B.  SYNOPTIC  CONDITIONS 

On  16  September,  a  sea  level  high  pressure  center  of  1031  mb  was  located  in  the 
eastern  North  Pacific  Ocean,  northwest  of  California,  along  with  a  weaker  high  prepare 
over  Arizona.  A  thermal  trough  was  present  in  the  Central  Valley  of  California  and  in- 
tensified during  the  day.  which  led  to  alongshore  flow  in  the  afternoon.  Winds  were 
generally  from  the  northwest  along  the  coast,  with  a  reported  wind  speed  of  5  m  s  from 
due  west  at  Monterey  in  the  morning.  The  wind  shifted  to  the  northwesterly  in  the  aft- 
ernoon. Skies  were  generally  clear  on  the  the  16th.  The  temperature  was  56°  F  at  A  a.m. 
local  and  rose  to  a  maximum  of  76°  F  later  in  the  day.  Fig.  1 1  on  page  38  shows  the 
synoptic  conditions  for  this  day. 

C.  WIND  FIELD  DATA  AND  ANALYSIS 

1.     Wind  Direction  and  the  Onset  of  the  Sea  Breeze 

The  wind  direction  cross-section  for  10  September  was  checked  for  a  strong 
change  in  wind  direction  with  time  that  indicated  the  passage  of  the  sea  breeze  front. 
The  onset  of  the  sea  breeze  occurred  at  1753  UTC  (0953  PST).  At  onset,  the  wind  di- 
rection veered  from  ISO0  to  300°.  The  temporal  wind  direction  shift  in  this  case  was 
90°  in  12  minutes.  The  wind  direction  cross-section  is  shown  in  Fig.  13  on  page  41. 
Following  this  directional  shift,  a  perturbation  in  wind  direction  occurred  directly  behind 
the  front.  The  wind  di  action  in  this  pertubation  at  50  meters  was  120°,  and  was  con- 
sidered to  be  the  offshore  flow  being  overriden  at  the  base  of  the  sea  breeze  front.  The 
presence  of  this  pertubaTion  may  possibly  be  explained  by  the  "lobe"  and  "cleft"  structure 
that  has  been  observed  in  the  three-dimensional  gravity  head  current.    This  pertubation 
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Fig.  11.      Synoptic  Situation  on  16  September:     Solid  lines  are  isobars. 

is  the  area  where  the  opposing  current  penetrates  the  gravity  current  head  and  causes  a 
"cleft"  in  frontal  boundary.    The  "lobe"  is  the  location  of  strongest  flow  at  the  front  of 
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the  longitudinal  vortex.  If  a  "cleft"  passed  the  sodar  site,  the  wind  direction  would 
change  rapidly  as  the  onshore  How  of  the  sea  breeze  was  interrupted  by  the  opposing 
flow  and  back  to  the  onshore  direction  as  the  sea  breeze  is  reestablished.  The  "lobe"  and 
"cleft"  structure  is  shown  in  Fig.   12. 

A  complete  directional  shift  of  larger  scale  occurred  about  16.78  decimal  day 
(1S43  UTC).  The  wind  direction  backed  with  decreasing  height,  and  reached  60°  at  50 
meters.  This  region  was  at  the  same  location  as  a  maximum  in  downward  vertical  ve- 
locity. A  possible  explanation  for  this  area  of  directional  shift  with  height  is  that  the 
descending  air  from  the  overturning  leg  of  the  longitudinal  vortex  affected  the  horizontal 
integrity  of  the  flow.  The  shift  in  wind  direction  occurred  within  the  confines  of  sea 
breeze  front  which  was  defined  by  the  directional  shear  region  at  its  boundaries.  The 
structure  of  the  vertical  velocity  field  which  complements  the  wind  shift  region  will  be 
discussed  later  in  section  4.3. 

The  height  of  the  directional  shear  zone  slowly  increased  with  time.   The  gradi- 
ent in  the  this  zone  was  180c  over  S3  meters  and  showed  that  the  layer  of  marine  air 
deepened  as  the  sea  breeze  became  well  established. 
2.     Wind  Speed 

In  Fig.  14  on  page  42.  the  increase  in  wind  speed  began  at  16.75  decimal  day 
(1800  UTC).  The  wind  speed  slowly  increased  from  2  m  s  to  3  m  s  in  about  30  minutes. 
The  offshore  flow  on  16  September  was  above  8  m  s  at  450  meters,  with  speeds  greater 
than  4ms  down  to  250  meters.  When  cross-sections  were  compared,  the  change  in  wind 
speed  occurred  immediately  after  the  change  in  wind  direction,  llelmis  (1987)  observed 
a  wind  speed  change  from  3  m  s  to  5  m/s  in  about  four  minutes  across  a  sea  breeze  front 
with  offshore  flow  of  5  m  s  and  a  change  from  2  m  s  to  3  m  s  over  15  minutes  with  flow 
parallel  to  the  shore. 

The  speed  shear  zone  dividing  the  onshore  and  offshore  flow  was  located  be- 
tween 250  and  350  meters.  Helmis  et  al.  (1987)  found  that  the  top  of  the  inflow  layer 
was  between  100  and  200  meters;  Simpson  et  al.  (1977)  observed  that  the  speed  shear 
zone  was  near  600  meters  at  the  head  of  the  sea  breeze  front  and  decreased  to  300  meters 
after  the  front  had  passed;  Lyons  and  Olsson  (1973)  observed  average  inflow  layer 
heights  of  500m,  ranging  between  100  and  1000  meters;  and  Ogawa  et.  al.  (19S6)  ob- 
served that  the  speed  shear  layer  was  between  150  and  320  meters. 

The  strong  offshore  flow  ( >  4  m  s)  was  above  250  meters  and  decreased  in  in- 
tensitv  with  time  as  the  height  of  the  sea  breeze  increased.    In  Fie.   15.  the  maximum 
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Fig.   12.       Three  Dimensional  Vie\>  of  a  Gravity  Current  Head:     The  "lobe"  and 
"cleft"  structure  is  shown,   (from  Simpson  et.  al.,1977) 

wind  speed  after  frontal  passage  was  10  m/s  at  500  meters  around  17.00  decimal  day 
(2400  UTC).  The  maximum  wind  speed  observed  by  Simpson  et  al.  (1977)  during  a  sea 
breeze  was  5  m  s  under  200  meters;  by  Johnson  and  O'Brien  (1973)  was  7.5  m  s  at  an 
elevation  of  300  meters;  by  Lyons  and  Olsson  (1973)  was  6.5  m  s  at  200  meters;  and  by 
Ogawa  et  al.  (19S6)  was  5.5  m  s  between  100  and  200  meters.  The  top  boundary  of  the 
sea  breeze  flow  was  not  as  clearly  evident  in  the  wind  speed  cross-section  as  it  was  in  the 
wind  direction  cross-section,  but  there  was  a  wind  speed  transition  around  550  meters. 
The  wind  maximum  of  8  m  s  above  this  level  was  associated  with  the  offshore  flow. 
3.     Vertical  Velocity 

In  Fig.  16  on  page  45,  a  maximum  upward  velocity  of  0.6  m/s  occurred  shortly 
after  16.75  decimal  day  (1800  UTC)  and  coincided  with  a  "bump"  in  the  wind  direction 
contour  plot.  This  "bump"  indicated  the  highest  vertical  extent  of  the  gravity  head 
structure  and  was  present  in  all  the  sea  breeze  fronts  that    passed    the    sodar    site. 
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LASBEX  WIND  DIRECTION  CROSS-SECTION 
SEPTEMBER  1987 
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Fig.   13.       Two-Hour  Cross-Section  of  Wind  Direction  (16  September  1987):     Plot 
begins  at  1730  LTC  through  1930  UTC. 

Maximum  upward  vertical  velocities  of  3  m  s  were  observed  by  Ogawa  et  al.  (19S6)  and 
of  1.5  m  s  by  Helmis  et  al.  (19S7).  The  top  of  the  vertical  velocity  maximum,  which  is 
marked  by  the  zero  vertical  velocity,  coincided  with  the  bottom  of  the  shear  zone  in  the 
direction  and  speed  cross-sections.  From  the  direction,  speed  and  vertical  velocity 
cross-sections,  the  elevation  at  the  head  of  the  sea  breeze  front  on  16  September  was 
between  ISO  and  220  meters.  The  depth  of  this  same  region  observed  by  Simpson  et. 
al.  (1977)  was  approximately  600  meters.    The  updraft  region  was  located  between  50 
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LASBEX  WIND  SPEED  CROSS-SECTION 
SEPTEMBER  1987 
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Fig.   14.      Two-Hour  Cross-Section  of  Wind  Speed  (16  September  1987) 

and  150  meters.    Helmis  et  al.  (1987)  observed  an  updraft  region  that  was  from  the  sur- 
face up  to  150  meters. 

A  region  of  downward  vertical  velocity  was  present  around  16.7917  decimal  day 
(1900  UTC).  The  maximum  downward  velocity  in  this  region  was  -0.4  m  s.  This  was 
consistent  with  descending  flow  in  the  rear  leg  of  the  longitudinal  vortex  within  a  gravity 
head  current.  Using  a  mean  advance  of  2.54  m/s  for  the  sea  breeze  fronts  during 
LASBEX,  the  horizontal  extent  of  the  longitudinal  vortex  on  16  September  was  8.1  km. 
This  dimension  is  similar  to  the  7  km  longitudinal  vortex  described  by  Simpson  et  al. 
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LASBEX  WIND  SPEED  CROSS  SECTION 
SEPTEMBER  1987 
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Fig.   15.      24-Hour  Cross  Section  of  Wind  Speed  ( 16  September  19S7) 

(1977).  This  vortex  is  an  oblong  circulation  that  is  located  fully  within  the  the  sea  breeze 
front.  The  vortex  observed  on  16  September  was  a  third  as  deep  as  the  vortex  observed 
by  Simpson.  This  difference  between  the  depths  of  the  vortices  was  probably  due  to  the 
strength  of  the  offshore  flow  on  16  September,  which  was  up  to  8  m  s.  This  strong  off- 
shore flow  flattened  the  front  and  restricted  the  depth  of  the  onshore  flow.  The  de- 
scending leg  of  the  longitudinal  vortex  observed  on  16  September  also  showed  up  in  the 
wind  direction  cross  section.  In  Fig.  13.  there  was  a  strong  temporal  gradient  of  the 
wind  direction  at  1S50  L'TC,  in  the  region  where  the  vertical  velocitv  has  a  cradient  of 
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0.4  m  s  over  5.6  minutes.    This  location  is  where  vertical  motion  changed  from  ascent 
to  descent. 

About  16.7917  decimal  day  (1900  LTC),  the  vertical  velocity  contours  indicated 
the  area  where  the  sea  breeze  was  reestablished  after  the  penetration  of  the  descending 
leg  of  the  longitudinal  vortex.  This  evidence  supports  Simpson's  theory  that  the  gravity- 
head  current  is  a  possible  description  of  the  actual  sea  breeze  frontal  structure. 
4.     Standard  Deviation  of  Vertical  Velocity 

The  location  of  the  strongest  fluctuations  in  the  vertical  velocity  is  shown  in 
Fig.  17  on  page  46.  The  highest  standard  deviations  of  vertical  velocity  occurred  at  the 
base  of  the  sea  breeze  front  and  around  the  upward  leg  of  the  longitudinal  vortex.  In 
fact,  the  whole  region  of  the  sea  breeze  front  showed  the  highest  values  of  standard  de- 
viation, indicating  strong  turbulence  and  vertical  velocity  fluctuations  in  this  area. 

The  behavior  of  a„  within  the  boundary  layer  has  been  found  to  follow  certain 
relationships  (Caughey  and  Palmer.  1979).  Since  the  o^  field  was  available,  this  re- 
lationship was  investigated.  In  Caughey  and  Palmer,  ct„  was  normalized  using  the 
scaling  velocity  w..  Scaling  velocity  in  the  convective  boundary  layer  is  determined  from 
the  heat  flux  based  on  the  following  equation  (Ogawa  et  al.,1986): 


w,  =  (4-^V/)1/3>  (4-1) 


where 

g  =  acceleration  due  to  gravity  (m.s-i). 

V  =  mean  potential  temperature  in  the  boundary  layer  (°C), 

w'6'0  =  sensible  heat  flu    at  surface  (///5_loC)and 

Zi=  depth  of  the  mixed  laver  (m). 

The  value  of  w,,  could  not  be  calculated  since  heat  flux  data  were  not  available.  How- 
ever, over  a  period  30  minutes,  the  change  in  heat  flux  was  considered  to  be  negligible 
and,  therefore,  the  value  of  w.  was  also  taken  to  be  constant.  The  value  of  z,  was  found 
as  the  height  on  the  aw  cross-section  which  was  indicative  of  the  base  of  the  inversion 
(i.e.  where  the  values  are  decreasing  quickly  towards  0). 
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Tiy.   16.       Two-Hour  Cross-Section  of  Vertical  Velocity  (16  September  1987) 

The  values  of  a,  were  plotted  versus  zjz,  to  compare  with  Ogawa  et  al.  (19S6). 
The  emphasis  of  this  comparison  was  on  the  distribution  of  o„.  In  Fig.  IS  on  page  47 
and  Fig.  19  on  page  4S,  the  maximum  value  in  ow  occurred  at  zjz,—  0.25  for  the  pre- 
frontal case  and  at  zjz  =  0.3  for  the  postfrontal.  The  heat  flux  was  allowed  to  van.',  so 
that  the  maximum  in  Caughey  and  Palmer's  (1979)  curve  matched  the  variance  of  the 
acoustic  sounder  values  of  am  are  at  a  maximum  in  comparison  with  zjz,  .  The  plot  by 
Ogawa  et  al.  (1986)  is  shown  in  Fig.  20  on  page  49. 
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LASBEX  SIGMA-W  CROSS-SECTION 
SEPTEMBER  1987 
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Fig.   17.       Two-Hour  Cross-Section  -  Standard  Deviation  of  w  (16  September  1987) 

Ogawa  observed  that  the  maximum  value  of  the  normalized  o„  occurred  at  z\z, 
=  0.3.  The  values  of  cw  decrease  with  height  as  Fig.  20,  which  contains  the  free  con- 
vection case  as  described  by  Caughey  and  Palmer  (1979).  Contrary  to  Ogawa,  the  values 
of  ow  in  the  LASBEX  plots  did  not  fall  off  as  dramatically  with  height  and  at  times  ac- 
tually increased  with  height.  This  was  because  of  the  strong  offshore  flow  above  400 
meters  on  16  September  caused  more  turbulence  above  the  inversion.  The  wind  profile 
that  Ogawa  observed  decreased  in  intensity  with  height  above  the  inversion,  which  gave 
lower  values  of  a... 
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Fig.  18.      Standard  Deviation  of  w  on  16  September:      ow  versus  z/z,  for  the  Pre- 
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5.     Horizontal  Dhergence 

The  divergence  cross-section  for  the  16  September  is  shown  in  Fig.  21  on  page 
51.  When  compared  with  the  wind  direction  cross-sections,  the  horizontal  divergence 
cross-section  showed  an  area  of  convergence  at  the  head  of  sea  breeze  front.  Fig.  21 
shows  that  the  maximum  values  of  convergence  occurred  at  the  head  of  and  within  the 
sea  breeze  front.  Values  of  convergence  in  these  regions  were  up  to  S  x  10~3s_1.  The  first 
region  is  located  at  the  top  of  the  maximum  in  vertical  velocity  and  the  second  is  located 
where  the  the  longitudinal  vortex  is  turning  back  into  the  foward  flow.  What  is  partic- 
ularly notable  is  that  this  cross-section  showed  that  the  convergence  followed  the  wind 
shift  boundary  of  the  sea  breeze  front  outlining  the  boundary  of  the  gravity  head  struc- 
ture, similar  to  the  divergence  Fields  of  Lyons  and  Olsson  (1973). 
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LASBEX  HORIZONTAL  DIVERGENCE  CROSS-SECTION 
SEPTEMBER  1987 
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Fig.  21.      Tnso-Houi-  Cross-Section  of  Horizontal  Divergence  (16  September  19S7) 

6.     Turbulent  Kinetic  Energy 

The  cross-section  of  the  TKE  for  16  September  showed  that  the  highest  values 
were  located  in  the  sea  breeze  front  where  mechanical  mixing  and  vertical  motion  gen- 
erated the  greatest  turbulence.  The  actual  maxima  in  turbulent  kinetic  energy  were  lo- 
cated in  the  region  where  the  overturning  air 'turned  back  into  the  forward  flow  and  at 
the  strong  shear  zone  near  the  base  of  the  front.  The  TKE  cross-section  for  this  day  is 
shown  in  Fig.  22  on  page  52. 
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Fig.  22.         Two-Hour  Cross-Section  of  Turbulent  Kinetic  Energy  (16  September 

19S7) 

D.     LIDAR  DATA 
1.     Color  Coding 

Color  coding  was  used  on  the  WPL  lidar  plots  to  show  the  intensity  and  direc- 
tion of  the  radial  velocity  field.  "Cool"  colors  were  used  for  velocities  towards  the  lidar, 
whereas  "warm"  colors  were  used  for  velocities  away  from  the  site.  The  colors  in  each 
case  were  graduated  with  the  "cool"  velocities  graduated  from  green  to  blue  to  purple. 
with    green    indicating    slower  speeds  and  purple  indicating  the  highest    speeds.    The 
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graduation  for  the  "warm"  colors  was  from  tan  to  yellow  to  red,  with  tan  being  the 
slowest  speeds  and  red  the  highest  speeds.  White  was  used  to  indicated  zero  radial  ve- 
locities. 

2.  RHI  Plots 

The  lidar  information  available  for  16  September  1987  closely  agreed  with  the 
wind  field  produced  by  the  sodar.  The  RHI  on  this  day  was  primarily  oriented  on  a  line 
of  90°.  The  separation  between  the  low  level  sea  breeze  and  the  intense  offshore  flow 
was  clearly  indicated.  The  depth  of  this  offshore  flow  started  at  500  meters  and  had  a 
depth  of  1200  meters.  The  flow  above  2  km  was  onshore.  About  1722  UTC,  there  was 
an  indication  that  the  sea  breeze  front  was  moving  away  from  the  lidar  site.  This  is 
shown  in  Fig.  23  as  an  increase  in  the  speed  of  the  flow  moving  away  from  the  lidar  to 
the  east  (right  side  of  picture)  at  the  lowest  levels. 

By  1735  L'TC,  the  sea  breeze  had  an  estimated  depth  of  about  200  meters,  based  on  the 
height  of  the  radial  velocity  field  in  the  lowest  layer.  This  agreed  with  the  depth  of  the 
onshore  flow  determined  by  the  sodar.  which  was  between  200  and  250  meters.  The 
actual  location  of  the  sea  breeze  front  was  not  determined  due  to  the  spatial  resolution 
of  the  output  information. 

There  was  one  RHI  view  of  the  sea  breeze  oriented  down  the  Salinas  Valley  of 
140&degree;  it  is  shown  in  Fig.  24.  The  onshore  flow  was  easier  to  detect  at  this  ori- 
entation because  it 

showed  up  as  a  deeper  and  stronger  flow  both  towards  and  away  from  the  lidar.  The 
speed  of  the  offshore  flow  from  the  Salinas  Valley  was  weaker,  with  a  radial  velocity  of 
4  m  s.  The  strong  flow  above  the  sea  breeze  ( >  6  m  s)  was  believed  to  be  due  to  the 
geography  of  the  area  and  is  explained  in  the  next  section. 

3.  PPI  Plots 

The  PPI  plots  give  a  better  indication  of  the  offshore  flow  region  and  the 
progress  of  the  sea  breeze  front.  The  5°  elevation  angle  plots  show  a  velocity  maximum 
at  520  meters  of  8  m/s.  This  is  shown  in  Fig.  25  as  a  velocity  maximum  towards  the 
lidar  from  the  northeast  and  the  corresponding  maximum  away  from  the  lidar  towards 
the  southwest.  The  election  and  orientation  of  this  radial  velocity  maximum  corre- 
sponded very  well  with  the  sodar  wind  direction  and  wind  speed  cross-sections.  A  region 
of  wind  speed  of  7  m/s  with  a  direction  between  60°  and  90°  was  present  above  the  sodar 
site  at  the  same  time  as  the  radial  velocity  maximum  was  above  the  lidar  site.  (Refer  to 
Fig.   13  on  page  41  for  the  wind  direction  and  Fig.   14  on  page  42  for  the  wind  speed.) 
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Fig.  23.      RHI  Lidar  Display  at  1722  UTC,  16  September  1987:     Top  picture  is  the 
RHI  display,  ["he  lidar  is  at  the  center.  Due  east  is  to  the  right  side  of  the 

picture.  Distances  anu  elevations  are  m  kilometers.  Velocities  are  m  in  s. 
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Fig.  24.       RHI  Lidar  Display  at  1805  UTC,  16  September  1987:     Top  picture  h 

RHI  display.  Down  valley  is  to  the  right  of  the  picture.  Distances  and  el- 
evations are  in  km.  Velocities  are  in  m  s 
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Fig.  25.  PPI  Lidar  Display  at  1716  LTC.  16  September  1987:  Elevation  Angle 
was  5°.  True  north  is  at  the  top  of  the  picture.  Rings  are  spaced  4  kilo- 
meters apart.  Velocities  are  in  m  s. 


55 


According  to  the  lidar,  the  offshore  flow  was  from  the  high  terrain  to  the 
northeast  and  from  the  gap  in  the  coastal  ranges  where  the  San  Juan  Valley  is  located. 
The  strength  of  this  offshore  flow  was  attributed  to  the  fact  that  the  San  Juan  Valley  is 
an  outlet  of  airflow  for  .he  higher  Santa  Clara  Valley  which  induced  a  strong  offshore 
current. 

The  PP1  plots  at  1°  elevation  showed  the  progression  of  the  sea  breeze  front. 
The  convergence  at  the  head  of  the  front  was  detected  as  it  moved  in  a  general  southeast 
direction  away  from  the  lidar  site.  The  lidar  showed  that  the  sea  breeze  front  was  in- 
fluenced by  the  surrounding  terrain.  At  1719,  the  sea  breeze  front  was  estimated  to  be 
located  at  or  around  the  lidar  site,  since  some  weak  convergence  was  approximately  1.5 
km  south.  As  shown  in  Fig.  26.  the  sea  breeze  flow  advanced  northeast  up  to  the  high 
terrain  approximately  7  km  due  east  of  the  lidar  and  north  to  the  mouth  of  Elkhorn 
Slough.  At  this  time,  the  sea  breeze  front  had  not  advanced  at  all  into  the  Salinas  Valley. 

The  high  terrain  to  the  east  effectively  blocked  the  further  advance  of  the  sea 
breeze  front  at  lower  levels.  The  calculated  height  of  the  terrain  at  7  km  was  about  120 
meters.  By  1743,  the  front  had  moved  about  4.S  km  south  and  southeast  from  the  lidar 
and  had  began  to  move  up  Elkhorn  Slough,  north  of  the  lidar.  This  is  shown  in 
Fig.  27  on  page  58. 

The  estimated  time  of  sea  breeze  frontal  passage  at  the  sodar  site  based  upon 
the  lidar  data  was  around  1800  LTC  and  compared  favorably  with  the  time  determined 
from  the  sodar  wind  held  data  .  which  was  1745  LTC.  Using  the  lidar  positions  of  the 
convergence  zone  at  the  sea  breeze  front,  the  calculated  speed  of  advance  was  2.09  m  s. 
This  compared  well  with  average  speed  of  2.54  m/s  calculated  from  surface  data. 
Fig.  2S  on  page  59  indicated  that  the  sea  breeze  had  completely  filled  the  mouth  of  the 
Salinas  Vallev  and  all  low  terrain  to  the  east  of  the  lidar  by  1S27  LTC. 
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Fig.  26.  PPI  Lidar  Display  for  1719  UTC,  16  September:  Elevation  Angle  was 
1°.  True  north  is  at  the  top  of  the  picture.  Rings  are  spaced  4  kilometers 
j  part.  Velocities  are  in  m  s. 
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Fig.  27.       PPI  Lidar  Display  for  1747  UTC,  16  September  1987:     Elevation  Angle 

was  1°.  True  north  is  at  the  top  of  the  picture.    Rings  are  spaced  4  kilo- 
meters apart.  Velocities  are  in  m  s. 
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Fig.  28.      PPI  Lidar  Display  for  1827  UTC,  16  September  1987:     Elevation  Angle 

was  1°.  True  north  is  at  the  top  of  the  picture.    Rings  arc  spaced  4  kilo- 
meters apart.  Velocities  are  in  m  s. 
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V.     SUMMARY  AND  CONCLUSIONS 

A.     MAPPING  THE  SEA  BREEZE  CIRCULATION 

1.  Accuracy  of  the  Sodar  and  Lidar 

Various  features  that  have  been  produced  in  laboratory  experiments  such  as  the 
flow  regions  of  a  gravity  head  current  were  observed  during  LASBEX,  most  notably  the 
structure  in  the  frontal  region  where  turbulence  and  mechanical  mixing  produce  large 
variations  in  velocity.  This  area  is  caused  by  the  overturning  of  the  longitudinal  vortex 
that  makes  up  the  major  How  within  the  frontal  region  and  its  interaction  with  the  am- 
bient Piows  surrounding  it. 

The  unique  combination  of  the  lidar  and  the  sodar  proved  to  be  an  effective 
method  to  detail  the  structure  of  the  sea  breeze  circulation.  The  instruments  both  per- 
formed well  in  a  marine  environment.  The  sodar  detected  and  mapped  the  passage  of 
the  sea  breeze  front,  whereas  the  lidar  showed  the  separation  between  the  onshore  and 
offshore  flow,  as  well  as  the  progress  of  the  convergence  zone  associated  with  the  sea 
breeze  front.' 

2.  General  Characteristics  of  the  Sea  Breeze  Circulation 

The  general  characteristics  of  the  sea  breezes  that  were  observed  during 
LASBEX  on  the  Monterey  Bay  Salinas  Valley  in  September  1987  were  similar  in  many 
ways  to  those  sea  breezes  seen  in  other  studies. 

The  wind  direction  change  over  time  was  found  to  be  an  accurate  indicator  of 
sea  breeze  onset.  The  times  of  onset  were  consistent,  occurring  around  1742  L'TC  at  the 
sodar  station,  which  was  6.5  km  inland.  Differences  in  the  relative  slope,  speed  and 
general  shape  of  the  sea  breeze  flow  was  observed  for  each  front  as  it  passed  the  sodar. 
The  directional  shear  in  the  vertical  meausred  by  the  sodar  gave  a  good  indication  of  the 
depth  of  the  onshore  flow.    The  lidar  presented  a  similar  depth  in  the  RHI  scan. 

The  observed  sea  breezes  advanced  onto  land  with  a  mean  speed  of  2.54  m  s. 
and  lasted  between  8  and  13  hours'.  The  maximum  height  that  the  onshore  flow  achieved 
averaged  about  659  met  rs  vertically.  The  slope  of  the  onshore  flow  measured  using  the 
height  of  the  directional  shear  zone  gave  values  that  ranged  between  1:5  and  1:17  in  the 
first  kilometer  horizontally  of  the  flow  and  between  1:56  and  1:220  in  the  established 
onshore  flow.  These  values  for  the  slope  of  the  onshore  flow  compare  favorably  with 
those  cited  by  Simpson  et  al.  (1977). 
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The  increase  of  wind  speed  and  the  maximum  updraft  region  both  occurred  after 
the  change  in  wind  direction.  This  agreed  well  with  the  location  of  the  longitudinal 
vortex  within  the  sea  breeze  front."  The  maximum  wind  speeds  of  the  sea  breeze  recorded 
at  the  sodar  were  from  6  to  10  m .  s.  Vertical  velocities  observed  during  the  sea  breeze 
in  the  mouth  of  the  Salinas  Valley  ranged  between  0.2  and  1.0  m/s  and  agreed  with 
typical  values  recorded  in  previous  studies,  such  as  Ogawa  et  al.  (19S6).  The  correlation 
between  depth  of  onshore  (low  and  the  vertical  extent  of  the  updraft  region  was  ob- 
served. A  larger  area  of  upward  flow  was  associated  with  weak  offshore  flow  and  a 
deeper  front.  The  location  of  downdraft  regions  corresponded  to  the  descending  leg  of 
the  longitudinal  vortex  and  gave  a  rough  measure  of  the  length  of  this  current.  On  16 
September,  the  length  of  such  a  flow  was  found  to  be  around  7  km  long,  and  compares 
well  with  the  vortex  observed  by  Simpson  et  al.  (1977). 

The  horizontal  divergence  calculations  made  using  the  vertical  velocity  field 
produced  values  of  convergence  that  at  times  were  up  to  12x10  35'.  These  were  found 
to  be  twice  as  large  as  the  values  calculated  by  Lyons  and  Olsson  (1973)  in  a  Chicago 
lake  breeze  study.  The  convergence  values  calculated  from  the  LAS  BEX  vertical  veloc- 
ity field  are  believed  to  be  the  most  accurate  for  the  sea  breeze  front  to  date. 

?.  Tracking  of  the  Sea  Breeze  Front 
The  lidar  was  able  to  detect  and  track  the  convergence  zone  associated  with  the  leading 
edge  of  the  sea  breeze  front.  The  derived  mean  speed  of  advance  agreed  favorabh  with 
speeds  derived  through  surface  measurements  at  both  the  lidar  and  sodar,  as  well  with 
speeds  seen  in  other  studies.  '  The  depth  of  the  offshore  flow  detected  by  the  lidar  was 
over  twice  as  large  as  the  onshore  flow,  which  is  the  typical  dimension  seen  in  other 
studies.  However,  the  speeds  in  this  offshore  flow  were  either  the  same  of  larger  than  the 
onshore  flow.  This  observation  is  contrary  to  those  studies  which  saw  that  the  speed 
of  the  oilshore  flow  was  roughly  half  as  fast  as  the  onshore  flow.'  In  PPL  the  lidar 
showed  how  terrain  affected  the  movement  of  the  sea  breeze  front,  such  as  halting  the 
advance  of  the  sea  breeze  front  and  the  movement  of  the  onshore  flow  into  Llkhorn 
Slough. 

B.     RECOMMENDATIONS 

Since  many  of  the  features  of  the  sea  breeze  flow  have  been  observed  by  the  sodar, 
further  studies  should  concentrate  using  this  versatile  instrument  to  complete  the  dy- 
namic structure  of  the  sea  breeze  front,  with  concentration    on    the  variations  of  the 
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velocity  components.     These  variations  would  help  understand  the  dimensions  and 
structure  of  the  turbulence  within  the  sea  breeze  front. 

As  far  as  improvements  in  the  mapping  of  the  sea  breeze  structure,  there  was  an- 
other sodar  which  operated  concurrently  with  the  PSU  sodar,  but  whose  data  was  not 
available  at  the  time  of  this  study.  The  combination  of  this  study  and  the  data  from  the 
other  sodar  would  give  a  better  view  of  the  sea  breeze  structure  and  refine  the  actual 
spatial  dimensions  of  its  associated  flows. 
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APPENDIX  A.     TWO-HOUR  WIND  DIRECTION  CROSS  SECTIONS 
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Fig.  29.      Two-Hour  Wind  Direction  Cross  Section  -  17  September  1987 
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Fig.  30.      Two-Hour  Wind  Direction  Cross  Section  -  18  September  1987 
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Fig.  31.      T»o-Hour  Wind  Direction  Cross  Section  -  19  September  1987 
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Fig.  32.      Two-Hour  Wind  Direction  Cross  Section  -  20  September  1987 
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Fig.  33.      Two-Hour  Wind  Direction  Cross  Section  -  21  September  1987 
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Fig.  34.      Two-Hour  Wind  Direction  Cross  Section  -  22  September  1987 
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Fig.  35.      T^o-Hour  Wind  Direction  Cross  Section  -  23  September  1987 
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Fig.  36.      T\\o-Hour  Wind  Direction  Cross  Section  -  27  September  1987 
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Fig.  37.      T^o-Hour  Wind  Direction  Cross  Section  -  29  September  1987 
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APPENDIX  B.     TWO-HOUR  WIND  SPEED  CROSS  SECTIONS 
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Fig.  38.      Tno-Hour  Wind  Speed  Cross  Section  -  17  September  1987 
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Fig.  39.      Two-Hour  Wind  Speed  Cross  Section  -  18  September  19S7 
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Fig.  42.      Two-Hour  Wind  Speed  Cross  Section  -  21  September  1987 
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Fig.  43.      Tn o-Hour  Wind  Speed  Cross  Section  -  22  September  1987 
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Fig.  45.      Two-Hour  Wind  Speed  Cross  Section  -  27  September  1987 
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APPENDIX  C.     TWO-HOUR  VERTICAL  VELOCITY  CROSS  SECTIONS 


83 


O 

LJ 
CO 

CO 

uo 
o 

(CJ 

o  r\ 

CO 
>H  en 


cj  ct 

O  LJ 
_J  CD 
LJ  2Z 

>  LJ 

_J  C_ 

d  LJ 
O  CO 


01 

LJ 


X 

LJ 
CO 
CO 


)ss    c-oos    cross-    o'oor    cos:    o-oc:    o'oss    o-oo?    o-Ds;    croc:  cros 
(Sd3I3W)    1H3I3H 


Fig.  47.      ~h\o-Hour  Vertical  Velocity  Cross  Section  -  17  September  1987 
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Fig.  -48.      Two-Hour  Vertical  Velocity  Cross  Section  -  18  September  1987 
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Fig.  49.      Two-Hour  Vertical  Velocity  Cross  Section  -  19  September  1987 
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Fig.  50.      Two-Hour  Vertical  Velocity  Cross  Section  -  20  September  1987 
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Fig.  51.      T»o-Hour  Vertical  Velocity  Cross  Section  -  21  September  1987 
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Fig.  52.      Two-Hour  Vertical  Velocity  Cross  Section  -  11  September  1987 
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Fig.  53.      Two-Hour  Vertical  Velocity  Cross  Section  -  23  September  1987 
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Fig.  54.      Tuo-Hour  Vertical  Velocity  Cross  Section  -  27  September  1987 
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Fig.  55.      Tno-Hour  Vertical  Velocity  Cross  Section  -  29  September  1987 
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APPENDIX  D.     TWO-HOUR  STANDARD  DEVIATION  CROSS 
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Fig.  56.      Two-Hour  Standard  Deviation  Cross  Section  -  17  September  1987 
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Fig.  57.      Tm o-Hour  Standard  Delation  Cross  Section  -  18  September  1987 
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Fig.  58.      Tuo-Hour  Standard  Deviation  Cross  Section  -  19  September  1987 
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Fig.  59.      Tno-Hour  Standard  Deviation  Cross  Section  -  20  September  1987 
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Fig.  60.      Two-Hour  Standard  Deviation  Cross  Section  -  21  September  1987 
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Fig.  61.      Two-Hour  Standard  Deviation  Cross  Section  -  22  September  1987 
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Fig.  62.      Tno-Hour  Standard  Deviation  Cross  Section  -  23  September  19S7 
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Fig.  63.      T«o-Hour  Standard  Deviation  Cross  Section  -  27  September  1987 


101 


CJ 
LJ 

CO 

I    r\ 

en  co 
CO  co 
o  — < 

CJ  cz 

LJ 


CJ  LJ 

CD  Q_ 
— ■  LJ 
CO  CO 

X 

LJ 
CD 
CO 

en 


'£?■■■ 


0-055    coos    crcsv    o'oos-    o'osl   o-oc:    cross    croo:    o'osi    o-oo;  o-de 
(Sa3I3U)    1HSI3H 


Fig.  64.      Two-Hour  Standard  Deviation  Cross  Section  -  29  September  1987 
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APPENDIX  E.     24-HOUR  WIND  DIRECTION  CROSS  SECTIONS 

A.     NOTES 

The  cross-section  for  23  September  1987  was  too  noisy  between  23.958  and  24.208 
decimal  daw   Data  in  this  region  were  considered  unusable. 
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Fig.  65.      24-Hour  Wind  Direction  Cross  Section  -  17  September  1987 


104 


HI 
«  HI 

■       cc 


CO 

CO 

CO 

On 
CC  oo 
O  ^ 

oS 

F=cq 
OS 

LU  LU 
fTr- 

Qgj 
QCO 


X 

LU 
CQ 
GO 
< 


0051     0099     0  095     0  09»     005C     009Z     0051    0'09 

(SH313W)  1HOI3H 


Fig.  66.      24-Hour  Wind  Direction  Cross  Section  -  18  September  1987 
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Fig.  67.      24-Hour  Wind  Direction  Cross  Section  -  19  September  1987 
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Fig.  68.      24-Hour  \Yind  Direction  Cross  Section  -  20  September  1987 
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Fig.  69.      24-Hour  Wind  Direction  Cross  Section  -  21  September  1987 
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Fig.  70.      24-Hour  Wind  Direction  Cross  Section  -  22  September  1987 
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Fig.  71.      24-Hour  Wind  Direction  Cross  Section  -  23  September  1987 
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Fig.  72.      24-Hour  Wind  Direction  Cross  Section  -  27  September  1987 
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Fig.  73.      24-Hour  Wind  Direction  Cross  Section  -  29  September  1987 
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APPENDIX  F.     24-HOUR  WIND  SPEED  CROSS  SECTIONS 

A.     NOTES 

The  cross-section  f>r  22  September  1987  was  too  noisy  between  23.958  and  24.083 
decimal  day,  where  the  speeds  are  greater  than  10  m/s.  Data  in  this  region  were  con- 
sidered unusable. 

The  cross-section  for  23  September  1987  was  too  noisy  between  23.958  and  24.208 
decimal  day.    Data  in  this  resion  were  considered  unusable. 
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Fig.  74.      2-4-Hour  Wind  Speed  Cross  Section  -  17  September  1987 
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Fig.  75.      24-Hour  Wind  Speed  Cross  Section  -  18  September  1987 
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Fig.  76.      24-Hour  Wind  Speed  Cross  Section  -  19  September  1987 
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Fig.  77.      24-Hour  Wind  Speed  Cross  Section  -  20  September  1987 
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Fig.  78.      24-Hour  Wind  Speed  Cross  Section  -  21  September  1987 
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Fig.  79.      24-Hour  Wind  Speed  Cross  Section  -  22  September  1987 
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Fig.  80.      24-Hour  Wind  Speed  Cross  Section  -  23  September  1987 
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Fig.  81.      24-Hour  Wind  Speed  Cross  Section  -  27  September  1987 
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Fig.  82.      24-Hour  Wind  Speed  Cross  Section  -  29  September  1987 
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Thesis 

c.l     The  sea  breeze  circula- 
tion during  the  Land  Sea 
Breeze  Experiment  (LASBEX) 
in  central  California. 


